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This thesis compares the current wholesale level consumable replenishment inven- 
tory model, now in use at Navy Inventory Control Points (ICPs), with a proposed Mean 
Supply Response Time (MSRT) Model. The purpose of the MSRT model is to intro- 
duce a readiness measure into wholesale inventory management. The objective of the 
MSRT model is to determine inventory depths which minimize the mean supply re- 
sponse time subject to not exceeding the inventory dollar investment provided by the 
CICP model for the same items. The MSRT model uses a marginal analysis optimiza- 
tion procedure. A comparative analysis of the models’ results indicates that the MSRT 
model provides consistently better supply system performance, in terms of supply mate- 
rial availabilitv (SMA) and mean supplv response time (in days), than the UICP model 
for items with a medium to low average quarterly demand. For medium to high quar- 


terly demand items, there is no significant difference in the models’ performances. 
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I. INTRODUCTION 


This chapter addresses the following issues: 
1. The motivation behind the research. 
2. The objectives of the research. 


3. A preview of subsequent chapters. 


A. MOTIVATION 

Since 1982 the Naval Postgraduate School (NPS) has been examining in depth the 
inventory control models in use at the Navy's Inventory Control Points (1CP). The two 
main ICP’s are the Ships Parts Control Center (SPCC) and the Aviation Supply Office 
(ASO). 

This work was undertaken at the request of the Naval Supply Systems Command 
(NAVSUP) in conjuction with NAVSUP’s Resolicitation Resvstemization Project. The 
primary objective of the resolicitation project was the acquisition of new computer 
hardware for the ICP’s, but the effort also provided the opportunity to make software 
improvements that would make better use of the new hardware. 

Initially, NAVSUP asked NPS to “develop improvements to the existing peactime 
Wholesale provisioning models for secondary items used by ... SPCC and ... ASO.” [Ref 
I: pp. 1] Provisioning models are those that deternune the range and depth of repairable 
and consumable items that should be initially procured in support of a new weapon 
system. A provisioning model's objective is to procure enough spares to provide ade- 
quate support for the Weapon svstem until the first wholesale replenishment order ar- 
rives. 

In response to NAVSLP'’s initial request, NPS proposed several provisioning models 
[Ref. 1] and conducted detailed comparisons between these models and those being used 
by the ICPs at that time. [Ref. 2] As a result of these analvses a new provisioning model 
was adopted at the ICP’s in 1984. This model is called the Mean Supply Response Time 
(MSRT) model because its objective is to minimize the MSRT for a set of secondary 
items belonging to a new weapon system. \fean supply response time 1s defined as the 
mean time it takes the supply system to satisfy the demand for an item. 

The logical next step was to develop a replenishment model having the same objec- 


tive function. This model would be used to make all subsequent buys of inventory after 


the initial provisioning buy. The initial effort in this direction concentrated on a re- 
pairable replenishment model. 

The initial attempt to extend the MSRT concept to repairable item management 
was described by Gormly in his master’s thesis at NPS. [Ref. 3] His thesis investigates 
an aggregate demand inventory model for repairable items that assumes: 

the probability distribution for the inventory position (defined as on-hand plus on- 
order plus in-repair minus backorders) was ... uniform with its equally likely states 
being a function of a weighted sum of the procurement quantity Qp and the repair 
batch induction quantity Qi. [Ref. 4: pp. 3] 
This average was then incorporated into the continuous review model of Hadley and 
Whitin, [Ref. 4: chap. 4]. A different approach was taken by Apple, also for his master’s 
thesis at NPS. [Ref. 5] The model presented by Apple took a queuing model approach 
towards evaluating the delay created by waiting for carcasses to accumulate for repair, 
and attritions (not repairable) to “accumulate” towards a procurement quantitv. The 
model was actually a provisioning rather than a replenishment model. 

In the spring of 1987, McMasters attempted to prove the inventory position prob- 

ability distribution assumed in Gormly’s thesis. He discovered that the probability dis- 


tribution assumption was incorrect. He also found that the P 


our 


repairable model 
developed by Fleet Material Support Office (F MSO) had been based in the same as- 
sumption and hence was also incorrect. MicMasters was able to develop a correct ver- 
sion of the probability distribution for the case of a repair survival rate of 100% (all 
carcasses inducted for repair will be successfully repaired). Research continues on the 
case of a repair survival rate less than 100°%. 

Fortunately, the problems with the probability distribution of the inventory position 
do not exist if only consumable items are considered. The inventory position distrib- 
ution is merely that of the exact formulation from Hadley and Whitin. [Ref. 4: chap 4] 
The consumable version of the MSRT model i:. also a limiting case of the repairable 
model. As a consequence, it is convenient to conduct a comparison between the con- 
sumable version of the MSRT model and the current UICP inventory model rather than 
waiting until the MSRT repairable model is completely developed. Such a comparison 
will provide insights into the benefits of the MSRT model and possibly result in its im- 


plementation for consumables rather than waiting for the repairable version. 


B. THESIS OBJECTIVES 


1. Develop a computer program for the current ICP model for the inventory man- 
agement of consumable items. 


2. Derive an MSRT model for consumable items and develop a computer program for 
it. 
3. Using ICP provided inventory management data, make a comparative evaluation 
of the performances of the proposed MSRT model and the current ICP model. 
C. PREVIEW OF SUBSEQUENT CHAPTERS 
In chapter II a discussion of the Navy Supply System is provided to give the reader 
an understanding of the environment in which the models are designed to operate. This 
is followed by an in-depth look at the current ICP model, including the mathematical 
formulas and the constraints that are applied. In chapter II] the MSRT model is pre- 
sented and the algorithim used for optimization is described. Chapter IV discusses the 
data, hardware and software used in the research. Chapter V presents the measures of 
effectiveness used to assess model performance, the results of the comparative analysis 
and an evaluation of these results. Chapter VI provides a brief summary and recomm- 


endations. 


H. THE NAVY SUPPLY SYSTEM AND THE UICP CONSUMABLE 
INVENTORY MODEL 


The intent of this chapter is to introduce the reader to the supply system in which 
the two inventory models that form the crux of this research are designed to operate. 
The chapter also explains how the UICP model functions within that system. This will 
be accomplished by: 


Giving an overview of the Navy Supply System. 
Describing the wholesale consumable inventory system. 
Explaining basic inventory theory. 


Deriving the equations used in the UICP model. 


SA ce a 


Discussing the constraints placed on the UICP model. 


A. NAVY SUPPLY SYSTEM 

The purpose of this section is to identifv the structure to which the Navy Supply 
System belongs and the components which make up the Navy Supply System. 

The Naval Supply Systems Command (NAVSUP) is responsible for the manage- 
ment of the Navy Supply System. The Commander of NAVSUP gives policy direction 
and guidance on inventory management of all material within the svstem. The actual 
wholesale supply system management functions are carried out by inventory control 
points (ICPs). The major Navy ICPs are the Aviation Supply Office (ASO) for aviation 
material and the Ships Parts Control Center (SPCC) for surface ship and submarine 
items, 

The ICPs provide two broad categories of support: program support and supplv 
support. While both categories are equally vital to the success of the Navy Supply 
System, the supply support category is of most importance to this thesis. Supply sup- 
port functions include: 

1. Budget development for parts support. 
Requirements determination. 
Material procurement. 
Integrated inventory management. 


Material distribution and issue. 


Dw bp wv 


. Repairables management. 


This thesis will concentrate on the functions of requirements determination and budget 
development for parts support. 

The Navy stock points are the holders of the inventories in the Navy Supply System. 
“The main mission of a stock point is the physical distribution of material.” [Ref. 6: pp. 


26] Their functions include: 


1. Receiving and stowing material. 


No 


. Issuing and shipping material. 


Ley) 


. Reporting receipts and issues to each item’s ICP. 


Stock points are located throughout the United States and overseas, near major con- 
centrations of operational and industrial customers. 

The two types of material managed by the Navv Supply System are repairables and 
consumables. “A repairable is an item of supply that can be made to function by a re- 
pair process after it breaks.” [Ref. 6: pp. 4] Examples of repairables include gear boxes 
and circuit boards. A consumable is an item which, when broken, during use, is not re- 
paired. Examples of consumables are valves, steam pipe and flanges. The inventory 
models involved in this thesis research deal only with consumables. 

The Navy groups similar items for administrative and management purposes into 
Cognizance Groups (COG), which are designated by two alphanumeric symbols or 
“digits”. The COG is used to identify the ICP or agency responsible for managing the 
item. It also, indicates whether the item is a repairable or consumable. The consumable 
items used in this research have a “COG” of 1H, indicating that the items are consum- 
ables managed by SPCC. 

To foster inventory management and flexibility, COGS are further divided by adding 
two more digits to them. These last two digits segregate items by essentiality or weapon 
system and by requisition frequency. For example, the 4-digit cog, IH4A represents a 
1H cog item having an essentiality of 4 (the highest) and a requisition frequency of 3 or 


more per quarter (denoted by the A). 


B. WHOLESALE CONSUMABLE INVENTORY SYSTEM 
Within the Navy Supply System there are three levels of inventory management of 


items: 


Wholesale - The inventory is managed by ICP inventory managers. The inventory 
managers have worldwide visibility and control of this inventory. 


Retail Intermediate - This level of inventory is designated to support a specific ge- 
ographic area and is normally managed by the Navy stock point in that 
geographic area. 


Retail Consumer - This is a level of inventory held for an activity’s own use. 


The inventory models in this research deal only with the wholesale inventory level. 
The general characteristics of wholesale level inventories are: 
1. Inventory levels are computed based on system-wide demand data. 


2. The ICP inventory manager knows where the material is located and has unre- 
stricted access to it. 


3. The ICP must approve all requisitions for items from these inventories. 
4. Material is assigned to one or more retail stock points for storage and issue by the 
ICP inventory manager. 

As stated in item 1 above, wholesale inventory levels are computed based on 
system-wide demand data. How does this demand data reach the ICP inventory man- 
ager? It is the responsibility of the stockpoints to submit the demand data to the ICPs. 
When a stock point issues a consumable to an end-use customer, the stock point reports 
the issue to the cognizant ICP via a transaction item report (TIR). Stock points also 
use TIRs to report receipt and redistribution of material. The information from the 
TIRs is inumediately processed into the ICP data files. This and other information in the 
files aids the ICP inventory manager in making decisions involving requirements deter- 
mination, material distribution and procurement of replenishment stock. Data from 
these files was used in the comparative evaluation of the current UICP and proposed 


MSRT inventory models. 


C. THE BASIS FOR THE ICP INVENTORY MODEL 
There are two major questions an ICP inventory manager must answer with regard 
to consumable items: 
1. How much to order; 1e., the reorder quantity (Q)? 


2. When to order; i.e., the reorder point (R)? 


In an unconstrained world the ICP inventory manager could order as much as the cus- 
tomer needed and order as frequently as demands occurred. Unfortunately, the real 
world is full of constraints with which the inventory manager must deal. These include 
inventory costs. 

The variable inventory costs which are dependent on Q and R are generally divided 


into three categories: 


Order Costs - costs accrued when determining and processing buys, and the 
stockpoints’ receipt and stowage costs. 


Holding Costs - includes costs of investment, storage, obsolescence and pilferage. 


Shortage Costs - costs of not having items when needed. aa 


These three costs plus the cost of buying the items (which is not dependent on Q and 
R) make up the total inventory costs. 

Two basic types of inventory models exist which consider these inventory costs when 
computing reorder quantity, Q and reorder point, R. They are the continuous review 
and periodic review models. The consumable model used by the Navy ICPs for whole- 
sale inventory is a tvpe of continuous review model called the min-max model. The 
min-max model seeks the values of Q and R which minimize the average annual variable 
inventory costs. 

In the min-max model, a replenishment order occurs when an item’s inventory po- 
sition (on-hand plus on-order minus backorders) hits or falls below R. At the time of 
ordering, the amount ordered :- Q plus the difference between the reorder point value 
and the inventory position value at that time. In the simplest form of the min-max 
model. Q is usually the amount which minimizes the average annual ordering and hold- 


ing costs. For consumables, this is the Wilson-Harris economic order quantity [Ref. 7]: 


§.1D 
OS te 
NéeIC 
where: 
A = administrative order cost; 
D = forecasted average quarterly demand; 


I = holding cost rate; and 

C = cost of one unit of the item. 

The reorder point, R, for the simplest form of the min-max model is determined from 
minimizing the costs of carrying safety stock and incurring backorders (shortages). It 


is a function of leadtime demand and the variability of demand. The formula for R is: 


R=(Dx L)+ SL; 


where: 


L = the forecasted mean procurement lead time in quarters; and 


SL = safety level, a function of demand and lead time variability. 
The reorder point is computed by first computing the optimal risk of a stockout 
according to the following formula [Ref. 7: pp. 145]: 


OIC 


Ses i) 


where: 


4 = shortage cost of one demand backordered for one year. 


The optimal risk is used with the probability distribution for demand during procure- 
ment leadtime to derive the actual reorder point value. This reorder point value is the 
mean leadtime demand (D x L which is designated by 4), plus the safety level. For the 


normal distribution, the safetv level is 
Sh = 7 


where: 
z = the standard normal deviate for the optimal risk; and 


o = the standard deviation of leadtime demand. 


How else is the performance of an inventory model judged besides by the minimi- 
zation of average annual variable costs? Typically, several other Measures of Effective- 
ness (MOEs) are used. A MOE is a function of the decision variables (in this case Q 
and R) and inventory system parameters. The two major MOEs the Navy Supply Sys- 
tem uses are: 


1. How well the available Navy Stock Fund (NSF) budget is allocated (another cost 
type MOE). 


2. Supply material availability (SMA)(customer service type of MOE). 


In this era of diminishing resource dollars, the Navy Supply System must allocate 
its NSF inventory investment dollars to maximize customer support. The UICP model 
does this by allocating these investment dollars to the optimal order quantities and re- 
order point which minimize total average annual inventory costs while attempting to 
achieve a minimum goal of 85% SMA. The dollar investment needed to support the 
inventory depth (order quantity + reorder point) is used by NAVSUP in formulating 


NSF budget requirements. 


D. THE UICP WHOLESALE CONSUMABLE INVENTORY MODEL 
In this section, the inventory model used by the ICPs to compute the reorder 

quantity, Q, and the reorder point, R, for consumables will be discussed. The following 
topics will be covered: 

1. Total variable cost equation. 

2. Model assumptions. 

3. Derivations of Q and R. 
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. Constraints. 


1. Total Variable Cost Equation 
The Department of Defense (DOD) provides policies for determining inventory 
levels in DOD Instruction 4140.39 and states that the objective is “to minimize the total em 
of variable order and holding costs subject to a constraint on time-weighted, 
essentiality-weighted requisitions short.” [Ref. 6: pp. 2] 
The total variable cost (TVC) equation in DOD Instruction 4140.39 looks like: 


TVC = Order Cost + Inventory Holding Cost + Shortage Cost ag 


The order cost, inventory holding cost and shortage cost are functions of Q and R. 


The LICP model. based on DOD Instruction 4140.39 is the minimization of an an- 
nual vaniable cost equation composed of the sum of three terms: 


Ordering Cost - an average number of orders per vear times the administrative cost 
of placing an order or contract. 


Holding Cost - average number of units of stock on hand at anv random point in 
time, times the cost to hold one unit of stock in inventory for a vear. 
—————— 


Shortage Cost - average number of requisitions on backorder at any random point 
in time, times the cost of not filling one requisition for a vear times the 
mulitarv essentialitv of the item. [Ref. 6: pp. 3-A-3] 


The total annual variable cost (TVC) equation is: 


n n 


N 
4D; Q; NEG 
TVC= QO. A; + 1¢(R, ap joa — Bj +BiQ;. R)) ra cc B{Q;. R)) 
\ t =m t 


=) i=] 


where: 


i = item index; 


= total number of items in a inventory; 


expected number of units of item 1 backordered at any random point in time; 


n 
B 
2 = shortage cost of one requisition backordered for one year, 
E = military essentiality of item i; and 

S 


= expected number of units per requisition for item 1. 


From reference 6 the formula for B(Q,, R,) is: 


(Qn R=, |" — RIE + Opa) — Fes ala 


where: 


F(+) = cumulative probability distribution of leadtime demand. 


2. Model Assumptions 
The ICPs use the following key simplifying assumptions in the LICP consum- 
ables inventory model. While these assumptions may not reflect the “real world” envi- 
ronment in which the inventory model operates, the assumptions are necessary to 
streamline and simplify the computation of Q and R. The assumptions are: 


1. Steady state demand environment exists (the demand distribution does not change 
over time). 


2. Infinite procurement quantities are possible. 
3. Continual review of the inventory position occurs. 


4. The quantity is sufficient to increase the inventorv position up to a value of 
O+R. 


Procurement cost, A, is a constant, independent of order quantity, Q. 
6. Demands and other model parameters associated with different items are inde- 
pendent. 
3. Derivations of Q and R 
a. Derivation of Order Quantity (Q) 
The derivation presented in this section was taken from reference 8. The 
optimal order quantity, Q, is determined by solving: 


éTvVc 44D , IC , CB(Q,R) JES 
SHEE 1D IC, BOM (30, HE) ag 
where: 
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Lod) 


B(Q.R) 
Q 


et [°c Os) Ia 


Lop) 


However, this 1s an equation involving an implicit function of Q and it is difficult to 
solve explicitly. It can be solved iteratively, but those calculations would be tedious for 


large inventory systems containing a large number of items. For this reason, the ICP’s 


ignore this term in 5 = 0 and solve for Q using the remaining terms. The result is 
§DA 
Cr09= = oa 


b. Derivation of Reorder Point (R) 


The optimal reorder point, R, is determined from 





emuC _ AE CBQ.R) \ _ 
aR sic+(ic+4 aa =) 


Limiting arguments to eliminate Q from this formula results in an optimal R satisfying 


SIG 


ee hk Sic + nee 


where: 


W = expected requisitions per quarter. 


4. Constraints 
The LICP consumable inventory model imposes several constrai' : on the op- 
timal Q and R from the total variable cost equations. These constraints are imposed 
because many “real world” factors cannot be included in the cost equations. These 
constraints will be discussed in this section. , 
a. Reorder Quantity Constraints 
The ICP’s impose constraints on the quantity they actually order so that 


their order quantity, Q, is determined from the following formula: 
QO = minimum{12D, maximum(Qeog, 1, D)}. 


This formula insures that: 


1. The reorder quantity is at least 1. 
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2. An item is reordered no more than once per quarter (by ordering at least D, the 
quarterly demand). 


3. The maximum size of an order is no more than 3 years’ demand (12D) per 
DODINST 4140.39. 
b. Reorder Point Constraints 

The reorder point is constrained in two steps. First, the nght hand side of 
the RISK formula is computed. That RISK value is then constrained to a value of no 
less than an ICP-set minimum risk and no greater than an ICP-set maximum risk. These 
risk limits are established to prevent under or over investment in any given item. 

The reorder point, R*, for the constrained risk is computed next. The dis- 
tribution used for anv item is dependent upon the demand characteristics of that item. 
CICP limits the choices of distributions to the normal, Poisson and negative binomial. 
The choice is based on a comparison of », where p» = Dx L, to an ICP-set parameter, 
called the probability breakpoint. If yw is less than the probability breakpoint a negative 
binomial probability distribution is used. If wis greater than or equal to the probability 
breakpoint a normal probability distribution is selected. A Poisson probability distrib- 
ution is used for very low-demand items, specially designated bv a “Mark Code of 0”. 

The reorder point, R*, is then further constrained to no less than zero or 
an ICP-set value, called the Numeric Stockage Objective (NSO is a minimum stockage 
‘evel established by an ICP on low demand items for insurance stockage purposes), no 
less than an ICP-set percentage of 4 or the number of wholesale stockpoints carrying the 
wholesale level for the item (called policy receivers) designated by the ICP, and no larger 
than an ICP-set percentage of # with on-hand assets not exceeding the shelf life quan- 
tity. The result is the final constrained reorder point, R : 

5. Selection of 4 
The shortage parameter, /, is chosen so that the resulting inventorv levels meet 
the SMA goal of 85%. Later it is reduced if sufficient funds are not made available (not 


a current problem under stock funding). 
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HI. MEAN SUPPLY RESPONSE TIME MODEL 


This chapter explains both the theory and the computational procedures used in the 


mean supply response time (MSRT) wholesale consumable inventory model. 


A. THE MSRT MODEL 

In the MSRT model the objective is to determine the depth of each item which 
minimizes the aggregate mean supply response time over all items within a given four- 
digit cog, subject to a budget constraint or maximum investment level for that four-digit 
cog. (As stated in chapter I, the mean supply response time is defined as the mean time 
it takes the supply system to satisfy the demand for an item.) 

Mathematically the MSRT model seeks the optimal maximum integer depths, 
Sie 01 = |, 2. ..., n)evhich minimize: 


S'D,x MSRT(SW}) 
MSRT=—>= 


n ’ 


subject to: 


n 
C; x SW, < InvesimeniLevel, 
— 


=i 


Where: 
1 = item index; 
D = forecasted average quarterly demand for item 1; 
SW = the depth of item i; and 


C = unit price of item i. 
It has been shown [Ref. 1] that: 


Time S,. 


MSRT;= De 


where TWUS stands for the expected time-weighted units short per year. It is defined 
[Ref. 4: pp. 185] as: 


SW, 
TWUS(SW)) = (un; — S Wi) + »y (SHV, — xi) x pkxp ud; 


xX, =1 


where p( x; » ) is the Poisson or negative binomial probability mass function for demand 
x during lead time and »= Dx L. If demand can be approximated by the normal dis- 
tribution then the value of p( x; uw ) for integer x values can be obtained by using a 
continuity correction on the normal. Therefore, the objective function can be rewritten 


as: 


A 
> TWUS{SW)) 
VERS = 


i=] 


and, because the denominator is a constant, the objective function reduces to: 


n 
> THUS(SIY). 

i=1 
Conceptually, time-weighted units short for an item are the number of demands that 
VE can't be satisfied bv issues from stock. weighted by the length of time that each demand 
remains unsatisfied (in a backorder position). In the case of consumables these demands 
normally remain backordered until a procurement arrives. 

The constraint in the MSRT model has been established, for comparison purposes, 
as the total investment level resulting from the UICP model for the same four-digit cog. 
That value is obtained by adding the Q, and R values to get the maximum UICP dzpth 
for each item, multiplying this depth by the unit cost C, , and summing these products 
Over alii — leet 


B. MSRT MODEL SOLUTION 
The optimization technique used on the MSRT model is marginal analysis because 
ae of the integer nature of demand and hence the depth SW. Basic to the marginal analysis 
approach is the ratio RR,, where: 


TOS CX — 1) — PCS (X 
Se LEN oo 
i 


This represents the ratio of the improvement in TWUS,, when the depth of item i is in- 
creased from X-1 to X, to the cost C; incurred. when the depth is increased by that one 
unit. 

To compute an optimal depth for the MSRT model using marginal analysis, all 
items are initialized with a depth of zero (SW=0) and RR is computed for X= 1 for all 
i. That item having the largest value of RR; is then increased by one in its depth (SW). 
Its RR value is next computed with X=2 and all RR, values are compared again. That 
item whose RR is largest has its SW increased by one unit, its X value is increased by 
one, and a new RR value is computed for it. After each unit increase in an item i, the 
left-hand side of the investment constraint is increased by C,. This procedure of com- 
paring RR, values, increasing one item’s depth and reducing the budget constraint is re- 
peated until the addition of the next unit would cause the cumulative investment level 
for the items within that four-digit cog to exceed the investment level the LICP model 
produced for the same four-digit cog. At this point the calculations are terminated. The 
calculations could be continued, however, for onlv those items whose unit prices remain 
low enough not to exceed the budget constraint on the next step. 

Each time an item goes through the procedure of computing the improvement 
brought about by adding one unit to its depth, its mean supply response time in days is 


also computed using the following formula: 


91x THUS(SIV) 


MSRT, = D 


If the item’s MSRT is less than 0.01 it is excluded from further depth increases. This 
value was set as a lower bound on anv item’s MSRT since further reductions would not 
significantly enhance system performance. The value of 0.01 is, however, arbitrary. 

As mentioned earlier, to produce an optimal depth, SW, using the MSRT model the 
depth of each item needs to be initialized at zero. If initialized at zero the model 
produced a depth of zero for many low demand items used in the comparative analyses. 
Although this is correct, depths of zero are rarely allowed in the LICP model. In the 
LICP model the O+ R depth is almost always one or greater because the order quantity 
is almost always constrained to be one and the reorder point is constrained to no less 


than zero. Despite the constraints in the UICP model for very slow movers this is not 
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an accurate representation of the way these items are actually managed. Usually these 
items have no material on-hand and no material on-order until a requisition is received. 
The true reorder point then is -1, the reorder quantity is | and the O+ R depth is Zero. 

In this research effort the MSRT model was run under three different conditions: 
depth initialized at zero; depth initialized at one; and depth initialized at the UICP 
computed reorder point. With the depth initialized at one the results may not be opti- 
mal, but they provide a basis for a comparison with the existing UICP model. 

The MSRT model was run with the depth initialized at the reorder point computed 
by the UICP model to assess the impact on the different measures of effectiveness and 
computer run times of never having a depths less than the UICP reorder points com- 
puted by UICP. 


IV. SOFTWARE, HARDWARE AND DATA 


This chapter describes the data that was used to run the models and the computer 


hardware and software utilized. 


A. SOFTWARE AND HARDWARE 

The original programming of the MSRT model was done in FORTRAN 66 by 
Professor Alan McMasters. The first step in this thesis research was to convert the or- 
iginal programming to WATFOR 77. The conversion was made to allow for easier 
program formulation, analysis (debugging) and structuring. The completed program is 
in WATFOR 77 Release 2.5. The operating system was VM CMS. The hardware on 
which the programming and running of the models was done was the IBM 3033 model 
mainframe at the W.R. Church Computer Center at the Naval Postgraduate School, 


Monterey, California. 


B. DATA 

The data that was used for the comparative analyses was from a tape containing the 
data for all the 1H cog National Item Identification Numbers (NIINs) on file at SPCC 
in the spring of 1985. The tape contained 113,647 1H records with each record pro- 
viding information on one NIIN. The records were in the CARES (Computation and 
Research Evaluation Svstem) format. 

Because of the program size and the limitation on virtual storage and CPU time 
available for this research, only select data was used. Computer runs were restricted to 
one hour of CPU time due to policies of the NPS Computer Center. Data samples were 
therefore selectively restricted in size to enable program processing within the one hour 
time limit. Data samples were also restricted in size because of the large number of ar- 
rays used in the program and the limited virtual storage available to an individual user. 

The process used for selecting data was to first compute the average quarterly de- 
mand for all items within each four-digit cog. Next, enough four-digit cogs were selected 
to represent the entire range of average quarterly demands. Finally, if all the NIINs 
within a given four-digit cog couldn't be processed through the models within time and 
space limitations, a subset of those items was selected which had an average demand 


closely approximating that of the four-digit cog’s entire population. 
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Table | is a breakdown of the entire 1H data set and information on the data sets 
that were used for the comparisons. The first entry in Table 1 shows the items from the 
data that were not assigned a four-digit cog. 


Table 1. BREAKDOWN BY 4-DIGIT COG OF THE DATA BASE 
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V. MEASURES OF EFFECTIVENESS 


This chapter discusses the Measures of Effectiveness (MOEs) used in the compar- 
ative evaluation and gives the computer results for each. The following is a list of MOEs 
used: 

1. Mean Supply Response Time (MSRT). 
. Supply Material Availability (SMA). 
3. Average Davs Delay (ADD). 
4. Average Davs Delay for Backordered Requisitions (ADDBO). 
5. Davs Of Safety Level (DOSL). 
6. CPL Run Time. 
7 benib Churn. 


to 


SMA (specifically SMA2), ADD and ADDBO are standard CARES (Computation and 
Research Evaluation System) MOEs [Ref. 9]. 


A. MEAN SUPPLY RESPONSE TIME (MSRT) 
Mean Supply Response Time is the mean time in days taken by the supply system 
to respond to the demand for an item. The aggregate MSRT for each four-digit cog is 


computed as follows: 


> [91 x THUS] 


MSRT=| — 
2, 


n 
i=] 


where: 
TWUS = Time-weighted units short in days for item i; 
D = Forecasted average quarterly demand for item i; and 


n = Number of items in the four-digit cog. 


Table 2 gives the results for the MSRT measurements for the MSRT and UICP models 
with the MSRT model depth initialized (ID) at zero, one and the LICP reorder point. 


Table 2. MODEL COMPARISONS FOR MSRT (IN DAYS) 
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I. MSRT Analysis and Observations 
As would be expected, given its objective function of minimizing mean supply 
response time, the MSRT model provides better mean supply response time than the 
UICP model in all cases. The improvement in performance becomes greater as the av- 
erage demand of the population decreases. 
The MSRT model with depth initialized at zero provides the best results, but the 
difference in performance between the three forms of the MSRT model is negligible until 


the average quarterly demand of the sample gets very small. 


B. SUPPLY MATERIAL AVAILABILITY (SMA) 

Supply Material Availability is the percent of all requisitions that are satisfied from 
on-hand stock over a year. The Chief of Naval Operations has set the current SMA goal 
at 85%. For purposes of this research SMA is computed two different ways, called 
SMAI and SMA2. SMA’ is a demand-based measure calculated as one minus the ex- 
pected number of backorders divided by the average quarterly demand. SMA2 is a 
requisition-based measure calculated as one minus the number of requisitions short per 
year divided by the average annual requisition frequency. The formulas for each are as 


follows: 
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Y (80; 
4x Sp, 
i=] 


SMAL=1- 


n 
> [RSPY] 
SMA2=1-| = —— 


4x) UW] 


i=] 


where: 
EBO = Expected number of backorders per year; 
RSPY = Requisitions short per vear; and 
W = Quarterly requisition frequency. 
Tables 3, 4, and 5 present the results for these two SMA MOEs for the MSRT 


model. with the depth for each item initialized at 0, 1 and the reorder point (RP), re- 


spectively. 


Table 3. MODEL COMPARISON FOR SMA WITH THE MSRT MODEL 
HAVING ID = 0 





Sa NIPPLE 
AVG. 


COG | DEMAND SMAI SMAI SMA2 SMA2 
99.32 
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Table 4. MODEL Ce FOR SMA WITH THE MSRT MODEL 
HAVING ID = 


SAMPLE 
AVG. UICP MSRT UICP MSRT 
COG | DEMAND SMAI SMAI SMA2 SMA2 














Table 5. MODEL ee FOR SMA WITH THE MSRT 
MODEL HAVING 
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AVG. UICP MSRT UICP MSRT 
COG | DEMAND SMAI SMAI SMA2 SMA2 
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1. SMA Analysis and Observations 
Relative to the UICP model, the MSRT model provides improved or equivalent 
SMA performance in all cases except for SMA2 in the case of very high demand items. 


For both measures of SMA the improvement in performance provided by the MSRT 
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model becomes greater as the average demand of the group of items processed becomes 
smaller. It is interesting to note that, in general, the performance of both the UICP and 
MSRT models deteriorates as the average demand of the population decreases. The 
better performance of the MSRT model then results from the fact that its SMA values 
decrease at a slower rate than those of the UICP model. 

Comparison of the performance of the three MSRT models, (depth initialized 
at 0, 1 or the reorder point) shows remarkably similar performance with very few differ- 
ences in SMA greater than +1 %. Again, the MSRT model with depth initialized at 
zero provides the best results. 


C. AVERAGE DAYS DELAY (ADD) 
Average Davs Delay measures how long on the average the supply svstem takes to 
fill a requisition. on the average aggregated over all requisitions. ADD for each four- 


digit cog is computed two different wavs as follows: 


n 
6 x Syren 
ADD| =——=——~— ; 
eS 
E x le 
n D, 
> x > {THEB] x ees i 
— Ht 
Ae —-— ; 
: x a 
1 
where: 


TWEB = Time-weighted average number of requisitions backordered per year for 
item 1. 


ADD1 is strictly a requisition-based calculation. ADD2 begins as requisition-based 


but is converted to demand-based by factoring in average requisition size. Tables 6 and 
7 give the ADD results. 
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Table 6. AVERAGE DAYS DELAY (ADD1) 


SAMPLE 
AVERAGE 
DEMAND 


1HOD 4567.95 27.69 27.79 meas 
iHoA | __si146 | 186 | 480 | 480 | 372 
236.03 





Table 7. AVERAGE DAYS DELAY (ADD2) 
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1. ADD Analysis and Observations 
The results of the ADDI MOE are slightly inconsistent. In the case of ADD1 
the MSRT model performs worse than the LICP model for high demand items and 
better than the LICP model for low demand items. Without exception, however, the 
MSRT models perform better than the UICP model for ADD2. Once again, the im- 


provement in performance increases as the average demand of the group decreases. The 
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performances of the three forms of the MSRT models are similar with the depth initial- 
ized at zero being the best. 


D. AVERAGE DAYS DELAY FOR BACKORDERED REQUISITIONS (ADDBO) 

Average Days Delay for Backordered Requisitions measures how long an average 
backordered requisition remains on backorder before it is filled. ADDBO is a standard 
CARES requisition-based MOE. The formula is: 


a0 = 


Table § gives the ADDBO results. 


Table 8. AVERAGE DAYS DELAY FOR BACKORDERED REQUISITIONS 
(ADDBO) 
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AVERAGE 
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736.03 














1. ADDBO Analvsis and Observations 
An quick analvsis of the results for this MOE demonstrates a clear advantage 
to the LICP model. 


ie DAYSIOPSAFETY LEVEL 

Days of safety level measures the amount of average stock on hand to protect 
against a stockout during leadtime. The higher the days of safety level the lower the risk 
of a stockout. The formula for days of safety level, which is computed by four-digit cog, 


1S: 


US 


365 x ) [Gx SL] 
DOSL=——>"———__ ; 


> [4x Gx DI 


f= 


where: 


DOSL = Cost weighted average safety level in days of supply at the forecasted mean 
quarterly demand rate of usage; 


SL = Safety Level (Item maximum depth — reorder quantity — leadtime demandin 
units for item i); and 


C = Unit price for item 1. 


THe UICP reorder quantity was also used for the reorder quantity of the MSRT model. 


Table 9 gives the days of safety level results. 


Table 9. DAYS OF SAFETY LEVEL (DOSL) 







SAMPLE 
AY ERAGE 
DEMAND 


1HOD 69.38 aT 
1HOA 511.46 138.07 146.81 146.81 138.37 
236.03 232.77 


















1. Days of Safety Level Analysis and Observations 
An analysis of the days of safety level MOE shows that the MSRT model with 
item depths initialized at either 0, 1 or the reorder point clearly outperforms the UICP 
model. As before, the degree of improvement increases as the average demand of the 


group decreases. 
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F. "€PU RUN TIME 

Central Processor Unit (CPU) run time measures how long the computer takes to 
compute the optimal depth for each model. CPU run time is an important MOE be- 
cause any additional run time may affect the ICP’s if their ADP capacity is already 


nearing its limits. Table 10 gives the CPU run time results. 


Table 10. CPU RUN TIME (CPURT) IN MINUTES 
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1. CPURT Analysis and Observations 

AS expected, the MSRT model with its marginal analysis procedure for deter- 
mining optimal item depths under a budget constraint requires significantly longer CPU 
run time than the LICP model. Initializing the depth of each item at the LICP com- 
puted reorder point decreased the MSRT model run time by about 50%. 

The run times encountered for the MSRT model during this research were sig- 
nificant. However. these can be reduced by use of VS FORTRAN, more efficient pro- 
gramming techniques and the more up-to-date hardware available at the ICPs. The run 
times are probably comparable to those encountered when a new 4 value is being 
computed selected by CARES for the LICP model. The value of 4 is not changed very 
often. The long runs to recompute maximum depths for the MSRT model need not 
occur very often either. Once the depths have been computed, each item’s MSRT value 
could be used and is available as a subsequent goal for fine tuning its depth on a quar- 


terly basis, if that is desired. 
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G. DEPTH CHURN 
To assess the changes in maximum depth that would be expected from implementing 
the model, the “depth churn” between the UICP model and the MSRT model with ID 


= 1 was computed. For item i, this churn is defined as: 


Churn, = UICP Model Maximum Depth, - MSRT Model Depth, 
=(Q,+ Ri) — SW; 


The churn data in Appendix I represents the four-digit cogs with the most extreme and 
least extreme churn. Each churn value is presented with its frequency of occurrence. 
Cogs 1HOA and 1HOD showed the most extreme churn while cogs IH3E and 1H4F 
showed the least. 
1. Depth Churn Analysis and Observations 

The characteristics of the items that showed the greatest degree of churn, one 
+ and one -—, for each of the four-digit cogs listed in Appendix I are presented in Tables 
11 and 12 below. The O value, which is the same for both models. is included along with 


each item’s unit price and quarterly demand. Depth is in units of inventory. 


Table 11. COGS WITH LEAST CHURN 


eco] She | A [ona |p [AE [SO 
COG DEPTH DEPT CHURN PRIGE DMD 
Garo | 2 | as 2 | tao0 | 1.61 







Keeping in mind that a negative churn indicates a greater depth for the MSRT 
model, it is clearlv seen that the MSRT model places greater emphasis on low-cost, 
high-demand items. This is the expected result given that the MSRT model criteria for 
increasing an item’s depth is based on a comparison of the cost-weighted ratios of im- 


provement in time-weighted units short. 


28 


Table 12. COGS WITH GREATEST CHURN 


UICP MSRT A UNIT QTRLY 
COG Bele a DEPTH CHURN Q PRICE DMD 


65083 


77 





VI. SUMMARY AND RECOMMENDATIONS 


A. SUMMARY 

This thesis is part of the Naval Postgraduate School's continuing effort to introduce 
more effective inventory models into the Naval Supply System’s management of whole- 
sale level inventories. The specific objective of this thesis was to develop an inventory 
model for consumable items having the objective function of minimizing the aggregate 
mean supply response time (MSRT) and then to compare the results of this model with 
the current consumable LICP model. The MSRT consumable model has the same 
readiness related objective as the recently implemented wholesale provisioning model for 
consumables. 

This objective was accomplished by developing the mathematics for the new model, 
developing computer programs for both models, processing ICP inventory management 
data through both models and conducting a comparative evaluation of the results using 
various measures of effectiveness. The main FORTRAN computer program is presented 


in Appendix A. The major subroutines are presented in Appendices B through H. 


B. CONCLUSIONS 
This research demonstrated conclusively that the MSRT model improves supply 
system performance for consumable items with no increase in investment. 
1. Standard Measures of Effectiveness Summary 
In addition to being superior relative to the aggregate MSRT measure of effectiveness, 
the MSRT model provides significant improvements over the UICP model in the fol- 
lowing supply svstem performance measures for items with a low, average quarterly de- 
mand: 
1. Supply material availability (SMA); 
2. Average days delay (ADD); and 


3. Days of safety level. 


The improvement provided by the MSRT model decreased as the average quarterly de- 
mand of the groups of items processed increased. Performance of the two models was 
equivalent for the items with a high average quarterly demand. This trend was consist- 
ently observed across all measures of effectiveness with the exception of Average Days 


Delav for Backorders. 
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2. CPU Run Time Summary 
The CPU run time required by the MSRT model with its marginal analysis op- 
timization procedure was significantly larger than the run times required by the LICP 
model. The amount of run time required for any group of items was found to be sensi- 
tive to the average quarterly demand of the four-digit cog’s population and the number 
of items in the cog. Since the focus of this research was an initial comparison of the 
performance of the two inventory models, programming efficiencies were not pursued. 
Reductions in CPU run time could be obtained by: 
1. Using VS FORTRAN instead of WATFOR 77. 
2. Using more modern computer hardware. 


3. Using more efficient programming techniques. 


C. RECOMMENDATIONS 

Replacement of the LICP consumable model with the MSRT model should be 
considered by the ICPs. While the long computer run times are a disadvantage of the 
MSRT model, it should not be n2cessarv to make such runs every quarter. Instead, each 
item’s depth can be adjusted to retain its same MSRT value as when the full four-digit 
cog optimization was last run. This approach is similar to that used to determine the / 
shortage cost. 

Viewing consumables as a limiting case of repairables, research should continue to- 
wards developing an MSRT model applicable to the wholesale inventory management 


of repairable items. 
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APPENDIX A. MAIN PROGRAM 


PROGRAM REPMOD 
TEST WHOLESALE CONSUMABLE REPLENISHMENT MODELS 

THIS PROGRAM RUNS BOTH THE UICP MODEL AND THE MSRT MODEL. 

THE UICP MODEL INITIALIZES X AT ZERO, ONE OR THE UICP COMPUTED 
REORDER POINT DEPENDING ON WHICH STEPS ARE COMMENTED OUT IN THE 
SUBPROGRAM "MODOPT" (STEPS 658 TO 670). THE PROGRAM ALSO COMPUTES 
THE CHURN BETWEEN THE TWO MODELS AND STORES IT IN A FILE THAT IS 
OPENED IN STEP 65. 

CHARACTER*23 NAME2/'MIN INVESTMENT LEVEL' / 

CHARACTER*23 NAME3/'MINIMUM MSRT '/ 

CHARACTER*23 NAME1/'UICP CONSUMABLES MODEL '/ 

CHARACTER*4 Q1/'1 '/,Q2/'4D '/ ,Q3/'E0Q "7,047 Do 
CHARACTER*4 RQ1/'1 '/,RO@2/"4D “/5RO03/' E00 4/)R05/. Dee 
CHARACTER*4 RQ5/'UICP'/,Q5/'UICP'/ 

CHARACTER*4 QM, QMR 

CHARACTER*1 PRIND( 1500) 

REAL BB(10),QBASIC,QTEMP 

INTEGER STOP( 1500) ,MARK( 1500), LOT( 1500) ,MD,PBP( 1500) , RLC( 1500) 
INTEGER N,NN,X( 1500) ,NSO( 1500) ,QMIN( 1500) ,NRPR( 1500) , XMSRT( 1500) 
INTEGER Y( 1500) ,QS,Q21,Q2C ,Q22,NPO,NOPT,KK,NI,ROP( 1500) 

INTEGER RESV( 1500) ,BO( 1500) ,PPRLT( 1500) ,PPR1( 1500) , PPR2( 1500) 
INTEGER PPR3( 1500) ,PPR4( 1500) ,AD( 1500) ,CD( 1500) ,OHRFI( 1500) 
INTEGER XCHURN( 1500) 

INTEGER OHNRFI(1500),PPR1YR( 1500) ,UPLIM, ROP1, XUICP( 1500) 

REAL D( 1500) ,PCLT( 1500), LAM( 1500) ,SL( 1500) 

REAL H,RISK,Z( 1500) ,ZN( 1500) ,T( 1500) ,Q( 1500) ,QQ( 1500) 

REAL C1(1500) ,HC( 1500) , POC( 1500) ,CHRNCST( 1500) , CHRNTOT 
CHARACTER*4 C0G1( 1500) ,COG2( 1500) 

REAL MODMST, MODNSF ,MSRTG( 10) ,MSRTGG ,BGT( 10) , PVAR( 1500) ,RST( 1500) 
REAL DMAD( 1500) ,RF( 1500) ,C11( 1500) ,E( 1500) ,RMIN( 1500) ,RMAX( 1500) 
REAL TOV( 2) ,QRI( 1500) ,QRR( 1500) ,QR( 1500) 

REAL QQR(1500),AS1( 1500) , SLC( 1500) ,OBSO( 1500) 

CHARACTER*4 COGG1( 1500) ,COGG2( 1500) 

REAL*8 B,BUDGG 

REAL RESLT,B19PT1,B19PT2,B19PT3,B19PT4, B19MIN, B19MAX 

COMMON SN(1500,9) 

EXTERNAL MODMST 

THE NEXT PARAMETER MUST BE SPECIFIED WHENEVER A NEW COG IS 
INTRODUCED. 

THIS NUMBER IS PROVIDED BY THE PRINTOUT FROM TEMPDATA PROGRAM. 
N=1500 

THE NEXT NUMBER SPECIFIES FULL TABLE LISTING(NPO=0) OR ONLY 

A SUMMARY (NPO=1) 

NPO=1 

NNN=0 

NN1=0 

THE VALUE OF Q MUST BE COMPUTED. THE PARAMETER MQ TELLS THE 
SUBROUTINE ORDQAN WNICH VALUE IS DESIRED; MQ=1 FOR Q=1, MQ=2 FOR 
Q=4*D, MQ=3 FOR Q=E0Q, MQ=4 FOR Q=D, MQ=5 FOR QUICP. 

QM IS EITHER Ql, Q2, Q3, Q4 OR Q5, RESPECTIVELY. 

MQ=5 


QM=Q5 
QMR=RQ5 
C ** NN IS THE NUMBER OF MEASURES OF 
C %%* PERFORMANCE; IF SMA AND MSRT ARE USED THEN NN=2. 
NOPT=2 
NN=2 
C *#* A COUNT OF VIABLE ITEMS MUST MADE (NOT HAVING STRANGE DATA). 
m0 
NI=0 
OPEN (UNIT=10, FILE='COA'" ,STATUS='OLD') 
OPEN( UNIT=11,FILE='CHRNOA’ ,STATUS="NEW' ) 
DO 23 I=1,N 
22 READ(10,898,END=24 )COG1(1) ,COG2(1) ,MARK(1),(SN(1,J),J=1,9),SLC(I) 
~PRIMO( 1) ,O8SO(1),E(1),PCLT(1),NRPR(I),C1(1),C11(1),DC1), 
*PVAR(I),RF(1),RESV(1I),BO(1),PPRLT(1),PPR1(1),PPR2(1),PPR3(1), 
*PPR4( I) ,OHRFI(1) ,OHNRFI(I),AD(I),CD(I), 
*AS1(1) ,DMAD(1),LOT(1) ,NSO(I) ,QMIN(I) 
S0¢mgeeat (2A2,11,9A1,1X,A1,1X%,A1,6X,F3. 2,F3. 3,F4. 2,7K,14,3F10. 2,10X, 
* Cron eno 5s , 17,1 ,17,1X,218,F8.0, 38%, F 10. 2,218,115) 
IF(D(1). GT. 250000. )GO TO 22 
TEGEC 1). LT.Ga0jGO TO 22 
IF(C1(1). GT. 999999. )GO TO 22 
iemeci¢ 1). LT:.01)GO TO 22 
IF(LOT(I).NE.0)GO TO 22 
IF(COG1(1).NE. '1H')GO TO 22 
KR-KK+1 
NI=NI+1 
23. CONTINUE 
24 DO 5 I=1,NI 
C0GG1(I)=C0G1(I1) 
COGG2( I)=C0G2(I) 
RMINCI)=. 10 
IF(COG2(1). EQ. '4A". OR. COG2(1). EQ. '3A'. OR. COG2(1). EQ. '2A'. OR. 
*CO0G2(1). EQ. '1A'. OR. COG2(1). EQ. 'OA' ) THEN 
RMAX(1)=. 30 
LAM(1)=1500 
PBP(1)=0 
RLC(I)=1 
ELSEIF(COG2(1I). EQ. '4B'. OR. COG2(1). EQ. '3B". OR. COG2(1). EQ. '2B'. OR. 
*CO0G2(1).EQ. '1B'. OR. COG2(1). EQ. 'OB' )THEN 
RMAX(1)=. 40 
LAM(1)=1000 
PBP(I)=0 
RLC(I)=1 
ELSEIF(COG2(1). EQ. '4D'. OR. COG2(1). EQ. '3D'. OR. COG2(1). EQ. '2D". OR. 
*C0G2(1). EQ. '1D'. OR. COG2(1). EQ. 'OD' ) THEN 
RMAX(1)=. 30 
LAM(1)=2000 
PBP(1)=0 
RLC(1I)=1 
ELSEIF(COG2(1). EQ. '4E'. OR. COG2(1). EQ. '3E". OR. COG2(1). EQ. '2E'. OR. 
*COG2(1). EQ. '1E'. OR. COG2(1). EQ. 'OE' )THEN 
RMAX(1)=. 40 
LAM(1)=2000 
PBP(1)=0 
REC CT )=! 
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ELSEIF(COG2(1). EQ. '4C'. OR. COG2(1). EQ. “30. OR COG2(1). EQs ze wmnHE 
RMAX(1I)=. 50 
LAM(I)=500 
PBP(1)=0 
REC i= 

ELSEIF(COG2(1I). EQ. '1C*. OR. COG2(1I). EQ. '0C' )THEN 
RMAX(I)=. 99 
LAM(I)=500 
PBP(1I)=20 
RLC(1I)=0 

ELSEIF(COG2(1I). EQ. '4F’. OR. COG2(I). EQ. '3F'. OR. COG2(1). EQ. '2F') THEN 
RMAX(1)=. 50 
LAM(1I)=750 
PBP(1)=0 
RLC(I)=1 

ELSEIF(COG2(1). EQ. '1F'. OR. COG2(1I). EQ. ‘OF' )THEN 
RMAX(I)=. 99 
LAM(I)=500 
PBP(I)=20 
Ric@hj=c 

ELSEIF(COG2(1). EQ. 'N1'. OR. COG2(1). EQ. 'N2'. OR. COG2(1). EQ. 'N3" )THEN 
RMAX(I)=. 35 
LAM(I)=1500 
PBP(I)=0 
REC(I)=1 

ELSE 
RMAX(1)=. 35 
LAM(1)=2500 
PBP( 1)=0 
REC(I)=1 

ENDIF 

STOP(1)=0 

IF(RF(I). EQ. 0. 0)RF(I)=D(1) 

IF(RF(1). GT. D(I))RF(1I)=D(I) 

IFCPCLIC!). GT. 18, )PCLIC1 )=18: 

IF( PCLT(1). EQ. 0. 0)PCLT(I)=5. 0 

IF(C11(1). GT. 999999. )C11(1)=999999. 

TECCUICT). LT, OCI T)=C ham 

IFCAS1( 1). LT. 0. OV ASIC 1 )=0. 0 

IF(AS1(1). GT. 250000. )AS1(1)=250000. 

IF(OBSO(1). LE. .01)OBSO(1I)=. 01 

IF(OBSO(1). GT. 1. O)OBSO(I)=1. 0 

MOC), gee 5 (Lob Ke) Cereal 

UPLIM= AMAX0( 100000, (NINT(16*D(1)))) 

PPR1YR( 1)=PPR1(1)+PPR2(1)+PPR3(1I)+PPR4(1) 

IF(OHRFI(1). LT. 0)OHRFI(1)=0 

IF(OHRFI(1). GT. UPLIM)OHRFI(I)=UPLIM 

IF(OHNRFI(I). LT. 0)OHNRFI(1)=0 

IF(OHNRFI(I). GT. UPLIM) OHNRFI( I )=UPLIM 

IF(AD(1I). LT. 0)AD(1)=0 

IF(AD(1I). GT. UPLIM)AD( I )=UPLIM 

IF(CD(1I). LT. 0)CD(1)=0 

IF(CD(1I). GT. UPLIM)CD( 1 )=UPLIM 

IF(BO(1). LT. 0)BO(1)=0 

IF(BO( I). GT. UPLIM)BO( 1)=UPLIM 

TP(RESV( 1) in, 0) RESVCL)—0 
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IF(RESV(I). GT. UPLIM)RESV(1)=UPLIM 
IF(PPRLT( 1). LT. 0)PPRLT(1)=0 
IF(PPRLT( 1). GT. UPLIM)PPRLT( 1)=UPLIM 
IF(PPR1YR(1). LT. 0)PPR1YR(1)=0 
IF( PPR1YR( 1). GT. UPLIM)PPR1YR(1)=UPLIM 
Z(1)=D(1)*PCLT(1) 
IF(Z(1). LE. 0. 0) THEN 
Z(1)=0. 0 
PVAR(1)=0. 0 
Goo 5 
ENDIF 
4 VIMR=150. 
VIMRP=PVAR(1I)/Z(I) 
IF(VIMRP. GE. VIMR)PVAR(1)=0. 0 
IF(PVAR(I).NE.0.0)GO TO 5 
PVAR(1)=4. 85%Z(1)**1. 5 
5 CONTINUE 
MSRTGG=0. 0 
CALL ICPMOD(NI,H,B,X,Z,C1,C11,D,QQ,NRPR,RF,PVAR, LOT, 
*DMAD ,E, LAM,MARK, PBP,SLC,RMIN, RMAX,RLC,NSO,AS1, 
*QQR ,QMIN,OBSO,PRIND,HC, POC) 
VSL=0. 0 
VAT=0. 0 
DO 10 I=1,NI 
ch Beis QQC 1) ;QQR( 1) 
SLCT)=X(1)-QQR(1I)-2(1) 
HUICP(1)=X(1) 
VSL=VSL+C11(1)*SL(I) 
VAT=VAT+C11(1)*D(1)*4. 
10 CONTINUE 
SLD=VSL*365. /VAT 
CALL PRTOUT(1,NAME1,QM,B,QQ,QRI ,QQR,N,NI,NN,X,Z,C11,D,NSRTGG, 
*COGG1,COGG2 ,NPO,TOV,QMR, PBP, PVAR, SLD, PCLT ,NNN,NN1) 
MOD=1 
CALL ICPSMA(NI,COGG1,D,RF,DMAD,X,QQ,QMIN, Z,PBP, PVAR,H, RMAX, LAM, 
*RMIN,MARK,SLC,PCLT,C1,E,MOD,QRI,NSO, 
*NRPR) 
MSRTGG=TOV( 2) 
BUDGG=B 
DO 41 I=1)NI 
IF(MQ. EQ. 5)Q(¢1)=QQ(T) 
IF(MQ. EQ.5)GO TO 41 
41 CONTINUE 
MOD=3 
CALL MODOPT(NI,NNN,B,MODMST,X,Z,D,QQR,C11,STOP,MSRTGG,MOD,PBP, 
*PVAR ,Q,NN1) 
wie COMPUTE ROP(I) 
DOl\g777 I= 1, NI 
ROP(I)=X(1)-Q(I) 
XMSRT(I)=X(I) 
7777 CONTINUE 
VSL=0. 0 
VAT=0. 0 
DO 43 I=1,NI 
SL(I)=X(1)-QC¢1)-2(T) 
IF(SL(1). LT. 0. 0)SL(1)=0. 0 
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VSL=VSL+C11(1)*SL(I) 
43 VAT=VAT+C11(1)*D(1)*4. 
SLD=VSL*365. /VAT 
CALL PRTOUT(3,NAME3,QM,B,Q,QRI,QQR,N,NI,NN,X,Z,C11,D,MSRTGG, 
*COGG1,COGG2,NPO, TOV, QMR, PBP, PVAR, SLD, PCLT, NNN, NN1) 
CALL ICPSMA(NI,COGG1,D,RF,DMAD,X,QQ,QMIN, Z, PBP , PVAR,H, RMAX, LAM, 
*RMIN, MARK, SLC, PCLT,C1,E,MOD,QRI,NSO, 
*NRPR) 
CHRNTOT=0. 0 
PRINT* 
WRITE (6,8184) 
8184 FORMAT('0','ITEM CHURN DMD UNIT PRICE CHURN PRICE') 
DO 8181 I=1, NI 
XCHURN(I) = XUICP(1I)-XMSRT(I) 
CHRNCST( 1)=XCHURN(1)*C11(1) 
CHRNTOT=CHRNTOT+CHRNCST( 1) 
WRITE(11,8186)XCHURN(1) 
8186 FORMAT (1X,18) 
WRITE(6 ,8183)1,XCHURN(I) ,D(1) ,C11(I) ,CHRNCST( I) 
8183 FORMAT('0',I4,18,F10. 2,F13. 2,F16. 2) 
8181 CONTINUE 
WRITE (6,8185)CHRNTOT 
8185 FORMAT('0',’SUM OF CHURN PRICES IS §$ ',F18. 2) 
*1500 CONTINUE 
99 STOP 
END 


*THE INPUT FOR THE READ STATEMENT OF THE MAIN PROGRAM IS IN THE CARES 
*(COMPUTATION AND RESEARCH EVALUATION SYSTEM) III FORMAT. 


i) 


*VARIABLE NAME DATA DEN(S ) CARES CC# 


*COG1 COG C003 1-2 
*COG2 CoG CO03W 3-4 
*MARK MARK CODE 5 

*SN STOCK NUMBER DO46D 6 - 14 
*SLC SHELF LIFE CODE C028 16 
*PRIND PROCUREMENT INDICATOR DO25E 18 
*OBSO OBSOLESCENCE RATE BO57 DSi= 27 
*E ESSENTIALITY CODE C008C 28 = 30 
*PCLT PROCUREMENT LEADTIME BO11A 31 - 34 
*NRPR NUMBER OF POLICY RECEIVERS a as) 
Cl REPLACEMENT PRICE B055 46 - 55 
Ci UNIT PRICE B053 56 - 65 
*2— QUARTERLY DEMAND BO74 66 - 75 
*PVAR PROCUREMENT VARIANCE BO19A 86 - 95 
*RF REQUISITION FREQUENCY A023B 96 - 105 
*RESV RESERVATIONS A013A 106 - 113 
*BO BACKORDERS AQ11 114 - 121 
*PPRLT PPRS DURING LEADTIME 122 - 126 
*PPR1 PPRS 1ST QTR AFTER LT 127 isi 
*PPR2 oN ae 132 - 136 
*PPR3 = 33RD 2 ake 137 = Wed 
*PPR4 oe OTH as > 142 - 146 
*OHRFI ON-HAND SYSTEM RFI 147 - 154 
*OHNRFI ON-HAND SYSTEM NRFI 155 - 162 
*AD AWARDED DUES 163 - 170 
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*CD 
*AS1 
*DMAD 
*LOT 
*NSO 
*QMIN 


COMMITTED DUES 

SET-UP COST B058 
DEMAND MAD A019(2)+A019A(2) 
DIPEROR TYPE QTY B070 
REORDER LVL LOW LIMIT QTY 8020 
MINIMUM PRODUCTION QTY B061B 
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171 
Ly9 
225 
235 
243 
ean 


178 
186 
234 
242 
250 
250 


APPENDIX B. SUBROUTINE ICPMOD 


C ‘ie THE CURRENT ICP CONSUMABLES MODEL 
C %% THIS SUBROUTINE COMPUTES ONLY THE PROCUREMENT Q AND ROP. 
SUBROUTINE ICPMOD(N,H,B,X,Z,C1,C11,D,Q,NRPR,RF,PVAR, LOT, 
*DMAD,E,LAM,MARK, PBP,SLC,RMIN,RMAX,RLC ,NSO,AS1,QQR,QMIN, 
*OBSO, PRIND,HC, POC) 
INTEGER N,NRPR(N) ,LOTCN) ,MARK(N) , ¥( 1500) ,R121,08,Q9, QMIN(N) 
INTEGER X(N) ,NSO(N) ,RLC(N) , PBP(N) ,ROP( 1500) ,ROP1,1,J,K 
CHARACTER*1 PRIND(N) 
REAL Z(N),C1(N) ,D(N) ,RF(N) ,PVAR(N) , DMAD(N) , AS1(N) 
REAL SL( 1500) ,RMIN(N) ,RMAX(N) ,C3( 1500) ,RISK( 1500) ,E(N) ,SLC(N) ,H 
REAL Q1,Q(N),LAM(N) ,C11(N) 


REAL QQR(N) ,AA,AC,RRCT,HC(N) ,OBSO(N) , POC(N) 
REAL TEMP1,TEMP2,T1,Q1A,Q1B,Q1C,Q1D,Q1E,Q1SQRT,Q1MAX 
REAL OLP1P2,TP1P2,21B23H, B19AOC, TSTRSK ,RISKTB( 50) , TVALU,RISKJ 
REAL RESLT, VIMTR, PVIMR, QRATIO, VRATIO, VLOG, PROB , SMPROB , B19PT3 
REAL B19PT1,B19PT2,B19PT4,B19MIN, B19MAX, B21MIN, B21PT1, B21PT2 
REAL*8 B 
B=0. 0 
weeteete RISKTABLE FOR THE NCRMAL DISTRIBUTION 
RISKTB(1) = 0.46017220 
RISKTB(2) = 0.42074060 
RISKTB(3) = 0. 38206860 
RISKTB(4) = 0. 33457830 
RISKTB(5) = 0. 30853750 
RISKTB(6) = 0. 27485310 
RISKTB(7) = 0. 24196370 
RISKTB(8) = 0. 21185540 
RISKTB(9) = 0. 18406010 
RISKTB(10) = 0. 15865530 
RISKTB(11) = 0. 13566610 
RISKTB(12) = 0. 11506970 
RISKTB(13) = 0. 09680050 
RISKTB(14) = 0.08075670 
RISKTB(15) = 0. 06680720 
RISKTB(16) = 0. 05479930 
RISKTB(17) = 0.04456550 
RISKTB(18) = 0.03593030 
RISKTB( 19) = 0.02871660 
RISKTB( 20) = 0.02275010 
RISKTB( 21) = 0.01786440 
RISKTB(22) = 0. 01390340 
RISKTB( 23) = 0.01072410 
RISKTB( 24) = 0.00819750 
RISKTB( 25) = 0.00620970 
RISKTB(26) = 0.00466120 
RISKTB( 27) = 0.00346700 
RISKTB(28) = 0.00255510 
RISKTB( 29) = 0.00186580 
RISKTB(30) = 0. 00134990 
RISKTB(31) = 0.00096760 
RISKTB(32) = 0. 00068710 
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eee COMPUTE HOLDING COST (CHC). 


RISKTB( 33) = 0. 00048340 
RISKTB(34) = 0. 00033690 
RISKTB(35) = 0. 00023260 
RISKTB(36) = 0. 00015910 
RISKTB(37) = 0. 00010780 
RISKTB(38) = 0. 00007230 
RISKTB(39) = 0.00004810 
RISKTB(40) = 0. 00003170 
RISKTB(41) = 0. 00002070 
RISKTB(42) = 0. 00001330 
RISKTB(43) = 0. 00000850 
RISKTB(44) = 0. 00000540 
RISKTB(45) = 0. 00000340 
RISKTB(46) = 0.00000211 
RISKTB(47) = 0. 00000130 
RISKTB(48) = 0. 00000079 
RISKTB(49) = 0. 00000048 
RISKTB(50) = 0. 00000029 
DO 6 I=1,N 


HC = STORAGE RATE(.01) + TIME 


AMRRKK PREFERENCE RATE(.10) + OBSOLESCENSE RATE (OBSO(I)) 
HC( I)=OBSO(1)+. 11 
xxxme COMPUIE AND CONSTRAIN THE RISK 
TEMP 1=HC(1)*D(1) 


TEMP2=RF(1I)*LAM(1I)*E(I) 


TRCREMP2 SEQ: 0. 0}) THEN 


T1=0. 


ELSE 


0 


T1=AMIN1(TEMP1*C1(1) /TEMP2, 999999. 0) 


ENDIE 


RISK(1)=T1/(T1+1. 0) 
RISK(I)=AMAX1(RISK(1) , RMIN(I)) 
RISK(1)=AMINICRISK(I) ,RMAX(I)) 


pe COMEUITE PROCUREMENT ORDER COST (POC) 


IF(MARK( I). EQ. 0. OR. MARK(I). EQ. 1. OR. MARK(1). EQ. 2) THEN 
POC(I)=660. 
ELSEIF(MARK(1). EQ. 3. OR. MARK( I). EQ. 4. AND. (8000**2*HC(I)* 


*( AMAX1((1-RISK(1)),(1-. 5)))/8*( 660. +AS1(1))).GT. (C1(1)*D(1)))THEN 


POC(I)=660. 


EGSe mac PRIND( 1)e@EQ. 0°. OR. PRIND(I).EQ. 'B’. OR. PRIND(1). EQ. '2') 


% THEN 
POC( 1)=1940. 
ELSE 
POC(1I)=1970. 
ENDIF 


wikis COMPUTE AND CONSTRAIN REORDER POINT (ROP) 
wisetes BASIC REORDER LEVEL *ivve% 


lowe 
aves 
aid 
woes 


dev 


THE TEMPORARY VARIABLE NAMES USED IN THE REMAINDER OF THIS 
SUBPROGRAM WERE USED IN ORDER TO MATCH AS CLOSELY AS POSSIBLE 
THE VARIABLE NAMES USED IN THE ACTUAL UICP PROGRAM 

OLP1P2=1. 0-RISK(1) 

TP1P2=RISK(I) 

Z1B23H=Z(1) 

B19A0C=PVAR( I) 

IF(MARK(I).EQ.0)GO TO 1500 


oo 


Seakdedede 
500 


600 


650 
700 


800 


Kicked 


1000 


1100 


1200 


1300 


KKK 


1500 
1510 


1520 


1540 
1550 


2500 


IF (Z1B23H. LT. PBP(I))GO TO 1000 


NORMAL DISTRIBUTION ‘oc 


TSTRSK=TP1P2 
IF(TP1P2. LE. 0.5)GO TO 600 
TSTRSK=1. 0-TP1P2 
DO 650 K = 1, 50 
IF(TSTRSK. LT. RISKTB(K))GO TO 650 
TVALU=1. 0*K 
IF (TSTRSK. EQ. RISKTB(K))GO TO 700 
J=K-1 
RISKJ=0. 5 
IF( J. NE. 0)RISKJ=RISKTB(J) 
TVALU=K-=( (TSTRSK-RISKTB(K) ) /(RISKJ-RISKTB(K) )) 
GO TO 700 
CONTINUE 
TVALU=50. 0 
IF( TP1P2. GT. 0. 5) TVALU=-TVALU 
TVALU=0. 1*TVALU 
K=0 
RESLT=Z1B23H 
IF(B19A0C. LE. 0.0)GO TO 2500 
IF(TP1P2. NE. 0. 5)RESLT=Z1B23H+( TVALU*SQRT( B19A0C) ) 
GO TO 2500 


NEGATIVE BINOMIAL DISTRIBUTION ***%% 


VIMR=1. 01 

IF(Z1B23H. NE. 0. 0) VIMR=AMAX1(VTMR, (B19A0C/Z1B23H) ) 
PVTMR=VIMR-1. 0 

QRATIO=PVTMR/VIMR 
VRATIO=Z1B23H/PVIMR 
VLOG=VRATIO*ALOG( VTMR) 

IF( VLOG. GT. 6.9)GO TO 500 

K=0 

PROB=EXP( -VLOG) 

SMPROB=PROB 

GO TO 1300 

K=K+1 

PROB=( ( VRATIO+( K-1)) /K)*QRATIO* PROB 
SMPROB=SMPROB+PROB 

IF( PROB. LE. 0. 0001)GO TO 500 
IF(SMPROB. LT. OLP1P2)GO TO 1200 
RESLT=1. 0*K 

GO TO 2500 


POISSON DISTRIBUTION ***#% 


IF (Z1B23H. GT. 6.9)GO TO 500 
K=0 

PROB=EXP( -Z1B23H) 
SMPROB=PROB 

GO TO 1540 

K=K+1 

PROB=( Z1B23H/K)*PROB 
SMPROB=SMPROB+PROB 

IF( PROB. LE. 0.0001)GO TO 500 
IF(SMPROB. LT. OLP1P2)GO TO 1520 
RESLT=1. 0*K 

GO TO 2500 

RESLT=RESLT 
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3860 


3900 
3920 


3940 
3990 


IF (D(I).GT.0.0)GO TO 3860 
RESLT=AMAX1( Z(1) ,FLOAT(NSO(I)),0. 0) 
CALL SHFLIF(SL(1I) ,SLC(I)) 

GO TO 3990 

CALL SHFLIF(SI(1) ,SLC(1)) 
B19PT3=SL(1I)*D(1)+Z(1) 

IF (NSO(I). GT. 100000. 0)NSO(I)=0. 0 
B19PT1=Z( 1)+33. *D(I) 
B19PT2=AMAX1(RESLT , FLOAT( NRPR(I))) 
TRCeSE( Is. LT. 7..0)GO TO 3900 
B19PT3=99999999. 0 

GO TO 3920 

B19PT3=B19PT3*D(1) 
B19MIN=AMIN1(B19PT1, B19PT2,B19PT3) 
B19PT4=1. 0*Z(1) 
B19MAX=AMAX1(B19MIN,B19PT4 , FLOAT(NSO(1)),0. 0) 
RESLT=B19MAX+0. 999 

ROP( I)=AINT(RESLT) 


6 CONTINUE 
wore COMPUTE Q NEXT 


4020 


*4100 
4100 


DO2 I=1,N 


Q1A=12. O*D(I) 
Q1B=D(1) 
Q1C=8. 0*( POC(1)+AS1(I))*D(I) 
QID=HC( 1)*C1(1)*AMAX1((1-RISK(1)),.5) 
Q1E=0. 0 
IF(Q1D. GT. 0. 0) THEN 
QIE=Q1C/Q1D 
ENDIF 
IF(Q1E. LE. 0. 0) THEN 
QISQRT=0. 0 
ELSE 
Q1SQRT=SQRT(Q1E) 
ENDIF 
QIMAX=AMAX1(Q1B,Q1SQRT) 
QI=AMIN1(QIMAX,Q1A) 
Q1=Q1+0. 999 
IF((SL(1)/4. ). LT. 7.0)GO TO 4020 
Q(1)=Q1+0. 999 
GO TO 4100 
B21PT1=AMAX1(ROP(1)-Z(1),0. 0) 
B21PT2=SL(1I)*D(1)-B21PT1 
B2 1MIN=AMIN1(B21PT2,Q1) 
Q(1)=B21MIN+0. 999 
IF(Q(1). LT. 1. 0)Q(1I)=1. 0 
QC I)=AINT(Q(I)) 
QQEC1)=Q01) 
X(1)=ROP(1)+IFIX(QQR(I)+0. 999) 
B=B+C11(1)*FLOAT(X(1)) 


12 CONTINUE 


ale 


RETURN 
END 


*VARIABLES NOT IDENTIFIED IN APPENDIX 1 


*VARIABLE NAME DATA 
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*N 


*B 
*X 


*LAM 
*PBP 
*RMIN 
*RMAX 
*RLC 
*QQR 
*HC 


NUMBER OF ITEMS BEING PROCESSED 

(NOT USED) 

BUDGET (INVESTMENT LEVEL) 

DEPTH 

ECONOMIC ORDER QUANTITY 

LAMBDA (GIVEN IN THE MAIN PROGRAM) 
PROBABILITY BREAK POINT ii 
MINIMUM RISK vd 
MAXIMUM RISK ay 
REORDER LEVEL CONSTRAINT 4 
ECONOMIC ORDER QUANTITY 

HOLDING COST 


APPENDIX C. SUBROUTINES/FUNCTIONS - MODOPT, MODMST 


AND GOALM 


* IN THIS APPENDIX THREE SUBPROGRAMS HAVE BEEN LISTED TOGETHER 
*BECAUSE THEY FUNCTION AS ONE SUBPROGRAM AND COULD HAVE BEEN WRITTEN 
*AS SUCH. 


C --- 


AQANIQANANAAN 
‘ 
’ 
‘ 


C --- 
C --- 


SUBROUTINE MODOPT(N,NNN,B,AMODEL,X,2ZN,D,QR,C11,STOP,GOALG,MOD,PBP, 


%*PVAR,Q,NN1) 


PERFORM OPTIMAL ALLOCATION FOR GIVEN MODEL USING 

MARGINAL ANALYSIS METHOD AND LOWER BOUNDING. 

N=NO. ITEMS 

B=INVESTMENT LEVEL OF STOCK FUND 

AMODEL=ENTRY POINT FOR MODEL TO USE ( STANDARDIZED ARGUMENTS ) 
=OPTIMAL ALLOCATIONS PER ITEM 

ZN= MEAN DEMAND DURING RESUPPLY TIME OR PPV 

C1=PROCUREMENT COST FOR EACH ITEM 

RR=WORK VECTOR TO STORE RATIOS 

INTEGER N,NNN,K,MK,STEP,X(N) ,STOP(N) , PBP(N) , INDEX( 1500) ,XL(1500) 

REAL ZN(N),QR(N) ,DCN) ,MR,RR( 1500) ,5R, TRY, PVARCN) ,QCN) 

REAL MODNSF ,MODMST,C11(N) ,GOALG 

REAL*8 B,BR,BUDLIM 

I=1 

SR=0. 

Kk=1 

NN1=0 

INITIALIZE 

BUDLIM=B 

B=0. © 

DO 2 I=1,N 

INITIALIZE STOP BEFORE OPTIMIZING ON INVESTMENT LEVEL(STOP=1 MEANS 

THAT THE LEVEL HAS HIT THE ITEM MSRT BOUND). 


WkAAKKTHESE FIVE STEPS PUT X AT ROP*#%%%% 


RP(1I) = X(1)-IFIXCQR(T)) 

fF CXC 1). GT. 1) THEN 
X(I) = RP(I) 

ENDIF 

B = B+(X(1)*C11(1)) 


scteiok to THESE POUR sohPs PUL X AT 1 OR OF Atakx 


we 
x 

te 

we 
oA 


* 


IF(X(1).NE. 0) THEN 
X(1)=1 
B=B+C11(1) 

ENDIF 


teexxethis step puts x at zero 


x(1)=0 

STOP(1I)=0 

INDEX(I)=0 

RR( I )=AMODEL( ZNC I) ,DC1) ,QRCI), ClIGL ex) +1,S1OP(1) .PBPC)), 
PVAR(1)) 


2 CONTINUE 
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12 
@ eas 
20 


30 


C --- 


C --- 
* 


* 


50 


@ sce 


10 


C--- 


We 


© sos 


STEP=0 
DO UNTIL MSRT OR BUDGET GOAL IS REACHED 
CONTINUE 
STEP=STEP+1 
MK=0 
MR=15000000. 
IF(MOD. EQ. 3)MR=0. 
DO 30 K=1,N 
IF(STOP(K). EQ. 1)GO TO 30 
IF( INDEX(I).EQ. 1)GO TO 30 
IF(RR(K).GE.MR)GO TO 30 
MR=RR(K) 
MK=K 
CONTINUE 
IF(MK .EQ. 0) GO TO 50 
ALLOCATE ONE MORE UNIT OF ITEM MK IF POSSIBLE. 
B=B+C11(MK) 
IF(B. GE. BUDLIM) THEN 
B=B-C11(MK) 
GO TO 50 
ENDIF 
X(MK)=X(MK)+1 
NEXT CHECK TO SEE IF GOAL HAS BEEN ATTAINED. 
CALL GOALM(X,N,ZN,D,QR,TRY ,GOALG, PBP, PVAR) 
IF(TRY. EQ. 0.0)GO TO 50 
SR=MR 
RR(MK)=AMODEL( ZN(MK) , D(MK) ,QR(MK) ,C11(MK) ,X(MK)+1, STOP(MK), 


*PBP(MK) , PVAR( MK) ) 


GO TO 20 
RETURN 
END 


ROUTINE TO MINIMIZE MEAN SUPPLY RESPONSE TIME -- MSRT 
REAL FUNCTION MODMST(ZZ,D,QR,C,K,STOP,PBP, PVAR) 
COMPUTE MARGIN ANALYSIS RATIOS ASSUMING 
POISSON DEMAND. 

REAL ZZ,QR,C,MSRT,MSRTD ,D,PVAR 

INTEGER K,STOP, PBP 

IF(ZZ.NE.0.0)GO TO 10 

MODMST=0. 0 

MSRT=-0. 0001 

ele) ae) a 

TS1=TWUS( Z2Z,QR,K-1, PBP, PVAR, MARK) 
TS2=TWUS( ZZ,QR,K, PBP, PVAR ,MARK) 
MODMST=(TS2-TS1)/C 

NOTE. MODMST WILL BE NEGATIVE 
MSRT=TWUS(ZZ,QR,K, PBP, PVAR,MARK) /D 

MSRTD=91. *MSRT 

IF(MSRTD. LT. 0. 01)STOP=1 

RETURN 

END 


ROUTINE TO SEE IF GOAL HAS BEEN ATTAINED. 
SUBROUTINE GOALM(X,N,Z,D,QR,TRY,MSRTG, PBP, PVAR) 
INTEGER N,X(N),XI,PBP(N) 
REAL Z(N),VO(1500),TV,SLT,D(N) ,QR(N) ,MSRTG, TRY 
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REAL MSRTC,MSRT, PVAR(N) 


TV=0. 0 
SLT=0. 
DO 10 I=1,N 
IF(Z(1). EQ. 0.0)GO TO 16 
XI=X(1) 
Str = ShueeeD(l) 
14 MSRTC=TWUS(Z(I),QR(I),XI,PBP(1) , PVAR(I) ,MARK)/D(T) 
MSRT=AMAX1(MSRTC, 0. 0) 
GO TO 18 
16 MSRT=0. 0 
18 VO(I) = 91. *MSRT 


TV = TV + VOCI)*DC(I) 


10 CONTINUE 


=TV/SLT 
TRY = 1.0 
IF(TV. LE. MSRTG)TRY = 0.0 


22 RETURN 
END 
* 


*VARIABLES NOT DEFINED IN PREVIOUS APPENDICES 
ve 


*VARIABLE NAME 
NNN 


DATA 

(NOT USED) 

MODMST 

LEADTIME DEMAND 

ECONOMIC ORDER QUANTITY 

1 OR O (INDICATES THAT THE ITEM HAS MET THE MSRT GOAL) 
MSRT GOAL (EQUALS THE MSRT THE UICP MODEL COMPUTED FOR 
THAT FOUR DIGIT COG 

(NOT USED) 

(NOT USED) 

LEADTIME DEMAND 

UNIT PRICE 

INTERMEDIATE DEPTH 

1 OR O (INDICATES THAT THE FOUR DIGIT COG HAS MET THE 

MSRT GOAL) 


45 


REAL*4 A/"A 
Us 
We 
i 
"7 ,A7/'7 
Sb. 
.A0)SL=100. 
.A0)GO TO 50 
.A)SL=1. /3. 
.A)GO TO 50 
mE Gla? 7 3. 
.B)GO TO 50 
16, UR, 216, 30). 
OO, SING, 180), 
.D)SL=4. /3. 
.D)GO TO 50 
.E)SL=5. /3. 
.E)GO 
ue: 


*E/'E 
*G/ 0G 
*N/'N 
*AL/'1 
*A6/"6 
REAL*4 


IF(SLC. 
IF(SLC. 
IF(SLC. 
IF(SLC. 
TFCSLG. 
TFCSLG; 
IF(SLC. 
IF(SLC. 
IF(SLC. 
IF(SLC. 
IF(SLC. 
IF(SLC. 
IF(SLC. 
IF(SLC. 
IF(SLC. 
IF(SLC. 
IF(SLC. 
TECSEG. 
IF(SLC. 
IF(SLC. 
IF(SLC. 
IFC(SLC. 
IF(SLC. 
IF(SLC. 
IF(SLC. 
IF(SLC. 
IF(SLC. 
IF(SLC. 
IF(SLC. 
IF(SLC. 
IF(SLC. 
IF(SLC. 
IF(SLC. 
IF(SLC. 
IF(SLC. 
IF(SLC. 
SL=100. 


50 RETURN 
END 


le 
ray 


*VARIABLE NOT DEFINED IN 


F 
G. 
nGs 
H 
H. 


Q 
Q. 
oR 
R 
S 
S 


APPENDIX D. SUBROUTINE SHFLIF 
SUBROUTINE SHFLIF(SL,SLC) 


'/,B/'B 


H/'H 
P/'P 
pO i 


SLC 


TO 50 


TO 50 
. SLC 
. oLc, 
. SLC. 
. SLC. 
7 oLG. 
. SLC. 


oe ye J 
oye 


We 


ESLC. EO: 
. SLC. EQ. 
PSLG. EG. 
SLC. EQ. 
/SIC.EQ, 
SLC. EQ. 
. J)SL=15. /3. 
.J)GO TO 50 
ok. OR) SILGO: 
UK. OR. SIG. EO 
Mh) SE=21 0/37 
»LYGO TG 50 
.M. OR. SLC. EQ. 
_M. OR. SLC. EQ. 
.N)SL=27. /3. 
.N)GO TO 50 
. P)SL=30. /3. 
. P)GO 


.45)GO TO 


EQ. 
.A7)GO TO 
. A8)SL=48. 
.A8)GO TO 
.A9)SL=60. 
.49)GO TO 


"/,C/'C 


¥A3/' 
A8 


A1)SL=3. /3. 
Al)GO TO 50 


A2)SL=6. /3. 
A2)GO TO 50 
A3)SL=9. /3. 
A3)GO TO 50 
A4)SL=12. 


A4)GO TO 50 


A5)SL=18. 
50 


A6)SL=24. 


A6)GO TO 50 


A7)SL=36 
50 


50 


50 
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1D) D9) nee 
K/L) eee 
R / XX J gay oes 
NTRS G7. 
AD'S "PRO 7 


Hse 


{3 


Hee 


es 


TE 
/3. 


PREVIOUS APPENDICES 


*VARIABLE NAME 
*SL 


DATA 
ITEM SHELF LIFE 
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APPENDIX E. REAL FUNCTION TWUS 


C --- ROUTINE TO CALCULATE THE EXP TIME WID UNITS SHORT FOR K UNITS 


10 


20 


21 


De, 


REAL FUNCTION TWUS(ZZ,QR,K,PBP, PVAR, MARK) 
REAL ZZ,P1,P2,SW,RP,QR,P3,P4,P5,P6 

REAL CCD1,CCD2,CCD3,CCD4 ,CCD5 ,CCD6, BETA1, BETA2 
REAL*8 PHI1,PHI2,DCD1,DCD2,D1,D2,T1,T2,Z,DBETA1,DBETA2,S,R,PV 
REAL*8 CD1,CD2,CD3,CD4,CD5 ,CD6 

INTEGER K,MARK,PBP 

SW=FLOAT(K) 

S=SW 

PV=PVAR 

a 

RP=SW-QR 

R=RP 

KRP=K-IFIX(QR+. 5) 

IF(ZZ. GE. FLOAT(PBP) )GO TO 20 

CALL CDFP(ZZ,KRP-1,P1,CD1,CCD1) 

CALL CDFP(ZZ,KRP,P2,CD2,CCD2) 

CALD CDEP (22,KRP+1.P3 0D3NCCD3) 

CALL CDFP(ZZ,K-1,P4,CD4,CCD4) 

CALL CDFP(2ZZ,K,P5,CD5,CCD5) 

CALL CDFP(2Z,K+1,P6,CD6,CCD6) 

IF (CCD1. LT. 0.000001) GO TO 10 
BETA1=(CCD1*ZZ%*2) /2. -CCD2%*ZZ**KRP+CCD3*KRP*(KRP+1) /2. 
IF(BETA1. LT. 0. 000001)BETA1=0. 0 

IF (CED4, GE. 0.000001) GO Ta 7 

BETA2=0. 0 

Corte 8 

BETA2=(CCD4*ZZ"*2)/2. -CCD5*ZZ*K+CCD6*K*(K+1) /2. 
IF(BETA2. LT. 0. 000001)BETA2=0. 0 

TWUS=( BETA1-BETA2)/QR 

RETURN 

TWUS=0. 0 

RETURN 

T1=(R-Z) /DSQRT( PV) 

TFL DABS( Ti), LE. 7. )GO TO 21 

PHI1=0. 0 

TFC iieGr 7. )Di=0, 6 

IF(T1. LT. (-7. 0) )D1=1. 0 

GO TO 22 

PHI1=(DEXP( -(T1**2)/2. ))/SQRT(2. *3. 14159265) 
CALL MDNORD(T1,CD1) 

D1=1. 0-CD1 

PRINT*, PV,D1,T1 

PRINT*, PHI1 
DBETA1=PV**(D1*%*(1. OFT 1°*2)-T1*PHI1)/2. 
T2=(S-Z) /DSQRT( PV) 

IF(DABS(T2). LE. 7. )GO TO 23 

PHI2=0. 0 

IF(T2. GT. 7. 0)D2=0. 0 

TF( 02 10-7 0) )D2=1),0 

GO TO 25 
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23 PHI2=(DEXP( -(T2**2)/2. ))/SQRT(2. *3. 14159265) 
CALL MDNORD(T2,CD2) 
D2=1. 0-CD2 

25 DBETA2=PV*(D2*( 1. 0+T2**2) -T2*PHI2) /2. 
TWUS=(DBETA1-DBETA2) /QR 
RETURN 
END 

* 


*ALL VARIABLES DEFINED IN PREVIOUS APPENDICES 
* 
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C wiv 


10 


il 


APPENDIX F. SUBROUTINE ICPSMA 


SUBROUTINE ICPSMA(N,COGG1,D,RF,DMAD,X,QQ,QMIN,Z,PBP, PVAR,H,RMAX, 
*LAM, RMIN,MARK,SLC,PCLT,C1,E, 
*MOD , QR,NSO,NRPR) 
INTEGER DEF,DEFR, ROP,ROPA,QA,QMIN(N) , XE, ZI ,Q2,X23,X24 
INTEGER XH,X2,X21,B019B,MARK(N) ,QS,Q21,Q2C,Y( 1500) ,X(N) ,MOD 
INTEGER NSO(N) ,NRPR(N) , PBP(N) , ROPR( 1500) 
CHARACTER*4 COGG1(N) 
REAL Z(N) ,DMAD(N) ,D(N) ,PVAR(N) ,PP( 1500) 
REAL TWEB( 1500) ,ETN( 1500) ,RF(N) , PCLT(N) 
REAL H 
REAL E(N) ,C1(N) ,SLC(N) ,QRCN) ,QQ(N) ,EF,EG 
REAL S1IZ( 1500) ,RMAX(N) , RMIN(N) , LAM(N) 
REAL*8 ZZ,PV,R,T1,PHI1,P1,CPHI1,USHRT1,WEB1,Q,T2,PHI2,P2,CPHI2 
REAL*8 USHRT2 ,WEB2,TR1,TR2,PHIR1,PHIR2,CPHIR1,CPHIR2,PR1,PR2 
REAL*8 EUS ,ERUS3,ERUS4, REPQ 
DO 100 I=1,N 
COMPUTE REQUISITIONS SHORT FOR PROCUREMENT CYCLE, ETN 
ZZ=2Z(1) 
PV=PVAR( I) 
ROP=X(1I)-IFIX(QQ(1)+0. 999) 
IF (D(1I).EQ. 0. 0) THEN 
ETN(I)=0. 0 
TWEB(I) = 0.0 
GO TO 100 
ENDIF 
IF(DMAD(I). NE. 0.0)GO TO 5 
DMAD(I)=(1. 39*D(1)**0. 75) 22 
GOMmens 
EF=1.57*DMAD(I)/(2. *D(1)) 
W1=0. 154 
EG=W1*D(I)*0. 5 
DEF=IFIX(EF+EG+0. 999) 
ROPA=ROP-DEF 
QA=IFIX(QQ(1)+0. 5)+DEF 
XE=MAX0( QA, QMIN(I)) 
Q=FLOAT( XE) 
IF(ROPA.GT.0)GO TO 10 
ZI=IFIX(Z(1)+0. 999) 
KN=MINO( ZI, XE) 
ETN(1)=(4. *RF(1)*FLOAT(KN) ) /FLOAT( XE) 
TWEB( I)=ETN(1)*(PCLT(1)+0. 5#W1)/4. 
GO TO 100 
IF(Z(1). LT. FLOAT(PBP(1)))GO TO 20 
R=FLOAT( ROPA) 
T1=(R-ZZ) /DSQRT( PV) 
Th( PARSON Gia7. colon 
PHI1=(DEXP( -(T1**2)/2. ))/SQRT( 2. *3. 14159265) 
CALL MDNORD(T1,P1) 
CPHI1=1.0 - Pl 
Go Te 12 
IF(T1. GT. 7. )CPHI1=0. 0 


50 


IF(T1. LT. (-7. })GPHI1=1. 0 
PHI1=0. 0 
12 USHRT1=(DSQRT( PV) )*(PHI 1-T1*CPHI1) 
WEB 1=PV*(CPHI1*( 1. +T1**2) -T1*PHI1) 
T2=(R+Q-ZZ) /DSQRT( PV) 
AB=(R+Q) 
IE@UASS(T2).GT. 7. )GO TO 13 
PHI 2=(DEXP( -(T2***2)/2. )) /SQRT(2. *3. 14159265) 
CALL MDNORD(T2,P2) 
CPHI2=1.0 - P2 
GO TO 14 
13 IF(T2.GT. 7. )CPHI2=0. 0 
IF(T2. LT. (-7. ))CPHI2=1. 0 
PHI2=0. 0 
14 USHRT2=(DSQRT( PV) )*( PHI2-T2*CPHI2) 
WEB2=PV*( CPHI2*( 1. +T2***2)-T2*PHT2) 
GO TO 50 
20 BOQ=PVAR(I)/Z(1) 
IF(BQ. GE. 1. 0. AND. BQ. LT. 1. 000001)BQ=1. 000001 
BK=Z(1)/(BQ-1. 0) 
PP( 1)=1. 0/( BQ***BK) 
USHRT1=Z( I) -FLOAT(ROPA)*(1. -PP(1)) 
WEB1=Z( 1) ***2+PVAR(I)-2. *Z( 1)*FLOAT(ROPA)+( (FLOAT(ROPA) )**2)*(1. - 
Peel) ) 
IF(ROPA. LE. 1)GO TO 23 
KM1=ROPA-1 
DO 22 J=1,KM1 
eet d+ )=( (BRI -1)/(J) )*CC80=1.))/BQ)*PP( J) 
USHRT1=USHRT1-FLOAT( J-ROPA)*PP( J+1) 
22 WEB1=WEB1-( FLOAT( J-ROPA)***2)*PP(J+1) 
23 KM2=ROPA+KE-1 
IF(KM2. LE. 1)GO TO 50 
USHRT2=Z( I) -FLOAT( KM2+1)**( 1. -PP(1)) 
WEB2=Z( 1 )**2+PVAR(1I)-2. *Z(1)*FLOAT(KM2+1)+( (FLOAT(KM2+1) )**2)*(1. - 
*EPC1)) 
DO 24 J=1,KM2 
Piel) —( (BKtJ-1)/(CI))*(( BOT. 1 /BQ)*PPCJ) 
USHRT2=USHRT2 -( FLOAT( J-KM2-1))*PP(J+1) 
24 WEB2=WEB2-(( FLOAT( J-KM2-1))**2)*PP(J+1) 
50 EUS=DMIN1((USHRT1-USHRT2) ,Q) 
ETN( 1)=EUS*4. *RF(I)/Q 
EW=(WEB1-WEB2)*RF( I) /D(1) 
TWEB( I)=EW/(2.*Q) 
100 CONTINUE 
C % COMPUTE AGGREGATE SMA 
101 SUMET=0. 0 
SUMTWB=0. 0 
SUMRF=0. 0 
DO 102 I=1,N 
SUNET=SUMET+ETN( 1) 
SUMTWB=SUNTWB+TWEB(I) 
102 SUMRF=SUMRF+RF(I) 
SMA2=1-(SUMET) /(4. *SUMRF) 
C we COMPUTE AVERAGE DAYS DELAY 
ADDA2=365. *( SUMTWB) /(4. *SUMRF) 
C * COMPUTE ADD FOR BACKORDERS 
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ADDBO2=ADDA2/( 1. -SMA2) 
SMA2=SMA2* 100. 0 


PRINT* 
WRITE(6 , 149) 
149 FORMAT('0',3X,'CARES REQUISITION-BASED PERFORMANCE MEASURES: ') 
WRITE( 6, 150)SMA2, ADDA2 , ADDBO2 
150 FORMAT('0” ,5X,'SMA=’ ,F8.4,5X, ADD= (f6uG,5%. ADDBO= Shion 
DO Wo 2-10 


IF(RF(1I). EQ. 0. 0)RF(1)=1. 0 
$1Z(1)=D(1) /RF(I) 
ETN( 1I)=ETN(1I)*SIZ(I) 
TWEB( 1I)=TWEB(1I)*SIZ(I) 

170 CONTINUE 
SUMET=0. 0 
SUMTWB=0. 0 
SUMRD=0. 0 
DO 180 I=1,N 
SUMET=SUMET+ETN( I) 
SUMTWB=SUMTWB+TWEB( 1) 

180 SUMRD=SUMRD+D( TI) 
SMA2=1-(SUMET) /(4. *SUMRD) 

C *%*%* COMPUTE AVERAGE DAYS DELAY 
ADDA2=365. *( SUMTWB)/( 4. *SUMRD) 
C %*%*% COMPUTE ADD FOR BACKORDERS 

ADDBO2=ADDA2/(1. -SMA”) 
PRINT* 
SMA2=SMA2*"100. 0 
WRITE( 6,190) 

190 FORMAT('0',3X,'DEMAND-BASED PERFORMANCE MEASURES: ' ) 
WRITE(6,150)SMA2 , ADDA2 , ADDBO2 
WRITE(6,194) 

194 FORMAT( ae ; Vo gedededeve de deve decleve de de vedesovere dese dese dese Fete se de Weede WIE Te sea Rese Ie eH ese ree | 
ae" Sevedevedese doves deve te clese vide de desde ae ae ae ae devotee ve deve vies a ae dese de ved desede dese dese sede ae ese ie dese ae deseo reacsess | ; 
ae * deresesevevesvee dese | ) 

RETURN 
* ; 


*VARIABLES NOT DEFINED IN PREVIOUS APPENDICES 
* 


*VARIABLE NAME DATA 

*COGG1 1H 

*QQ REORDER QUANTITY 
as 


END 
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APPENDIX G. SUBROUTINES - PRTOUT, OBJECT AND EBO 


*THE FOLLOWING THREE SUBPROGRAMS WERE LISTED TOGETHER BECAUSE 
*FUNCTION TOGETHER 


C --- ROUTINE TO PRINT OUT RESULTS 
SUBROUTINE PRTOUT(MD,NAME ,QM,B,Q,QRR,QR,NT,N,NN,X,Z,C11,D,MSRTG, 
*COG1,COG2,NPO,TOV,QMR,PBP,PVAR,SLD, PCLT,NNN,NN1) 
INTEGER ROP( 1500) ,NPO,NT,KQ( 1500) ,KQR( 1500) ,KQW( 1500) ,X( 1500) 
INTEGER PBP(N) 
REAL Z(N),0V(2,1500) , TOV(2) ,C11(N) ,MSRTG, BSW( 1500) ,SLD 
REAL Q(N),D(N) ,QRR(N) ,QR(N) , PVAR(N) 
REAL PCLT(N) 
REAL*8 B ‘ 
CHARACTER*23 NAME 
CHARACTER*4 COG1(N),COG2(N),QM,QMR 
COMMON SN(1500,9) 
DO 1 I=13N 
KQCI)=IFIX(Q(1)) 
ROP( I)=X(1)-KQ¢(1I) 
BSw(1) = €11(1)*X(1) 
1 CONTINUE 
CALL OBJECT(X,N,NN,Z,D,QR, PBP, PVAR, OV, TOV) 
WRITE(6 , 900) 
900 FORMAT( 1 x ; / / jf ‘ y fevevededere de tevede deve ete vevete te sevete te tevetee tele te dedete tet Teeter eee nete | 
ate! decedesededede dete de te tedete Fe eR IRE RET III IMI INN M ; 
ate V totactantaatantentatetoutaatertortertoats t 
WRITE(6 ,901)MD, NAME ,COG1(1) ,COG2(1) ,MSRTG,QM,N,NT,NNN 
OOImnORMATC O' 1X SMODEIC’ ,11, )',1X,A23,1X,'COG: " ,2A2,3x, 
* MSRP GOAL: ',FS,2,' DAYS’ ,3X, OP: ',A4,3X,'NI/N: ,14, 
ef 14.9%. LOWER: BD CODE: ,12) 
IF(NPO. EQ. 1)GO TO 907 
WRITE(6,902) 
902 FORMAT('0' ,5X,'NIIN',4X,'DEPTH' ,3X, 'MSRT(DAYS)' , 2X, 
Tames LVL. ,5X4 Gl1CS) ,2X,'Q *,3X, ROP’, 
covery. 8x, 'D' ,5X,' PCLT’ , 
*5X,'PVAR' ) 
WRITE(6,903)((SN(I,J),J=1,9),X(1) ,OV(2,1),BSW(I),C11(1),KQ(I), 
sROrPC 1). 2(1),D(1),PCLT(1),PVARC1),1=1,N) 
GOs erOnmi1( 3h 9A1,2X,15,2X,F10.4,1X%,F13. 0,1, F1i. 2,15,15, 
Silees Fs, 3,1X%,2X,F5.2,F9. 1) 
906 WRITE(6,904)TOV(2),B,TOV(1),SLD 
904 FORMAT('0',2X,'TOTAL PERFORMANCE: ' ,F9.3,1X,'$' ,F12.0,5X, 
*'SMA: ',F8.2,4X,'SAFETY LEVEL DAYS: ' ,F8. 2) 
GO TO 909 
907 WRITE(6,908)TOV(2),B,TOV(1),SLD 
908 FORMAT('0' ,4X, OVERALL PERFORMANCE: ',2X,F9.3,4X,'$',F12.2,5X, 
*'SMA: ',F8.2,4X,'SAFETY LEVEL DAYS: ' ,F8. 2) 
909 RETURN 
END 
% --- ROUTINE TO COMPUTE THE OBJECTIVE FUNCTIONS FOR GIVEN ALLOCATION 
SUBROUTINE OBJECT(X,N,NN,Z,D,QR, PBP, PVAR,OV, TOV) 
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* 


16 
18 


10 


10 


INTEGER N,NN,X(N),XI,PBP(N) 
REAL Z(N),OV(NN,N) ,TOV(NN) ,SLT,MSRT,OV1,D(N), QR(N) 
REAL CD,P,MSRTC,BO,BOT,DD, PVAR(N) 
TOV(1)=0. 
TOV(2)=0. 
BO=0. 0 
BOT=0. 0 
CONTINUE 
SLT=0. 
DO 10 I=1,N 
IF(Z(1). EQ. 0.0)GO TO 16 
XI=X(I) 
SLT = SLT + D(I) 
Ovi )=c: 
CALL EBO(2(1),XI,D(I),QR(1),0V1, PBP(I),PVAR(I)) 
MSRTC=TWUS(Z(I),QR(1),XI, PBP( 1) , PVAR( 1) ,MARK) /D(I) 
MSRT=AMAX1(MSRTC, 0. 0) 
OV(1,1)=AMIN1(OV1, 1. 0) 
GO TO 18 
OV(1,1)=0. 0 
MSRT=1. 0 
OV(2,I) = 91.*MSRT 
TOV(2) = TOV(2) + OV(2,I)*D(T) 
BO=B0+0V(1,1) 
CONTINUE 
IF(SLT. EQ. 0. 0)SLT=1. 0 
BOT=BO/SLT 
TOV(1)= (1. -BOT)*100. 
TOV(2)=TOV(2)/SLT 
RETURN 
END 


SUBROUTINE EBO(Z,X,D,QR,OV1, PBP, PVAR) 

REAL Z,D,QR,OV1,ALPHA1, ALPHA2 ,X1,X2,D1,D2,P1,P2,PVAR 
REAL*8 T1,T2,PHI1,PHI2,CDD1,CDD2,DD1,DD2,DLPHA1,DLPHA2,PV,CD1,CD2 
INTEGER PBP,X 

K1=X-IF1X(QR+. 05) 

X1=FLOAT(K1) 

X2=FLOAT(X) 

K2=X 

PV=PVAR 

IF(Z. GT. FLOAT(PBP))GO TO 10 

CALL GDEP(Z,K1,P1,GoinD1) 

CALL CDFP(Z,K2,P2,CD2,D2) 
ALPHA1=D1*(Z-X1)+X1*P1 

ALPHA2=D2*( Z-X2)+X2*P2 
IF(ALPHA1. LT. 0. 0)ALPHA1=0. 0 
IF(ALPHA2. LT. 0. 0)ALPHA2=0. 0 
OV1=(ALPHA1-ALPHA2)*D/1 

ele) dio) 20 

T1=(X1-Z) /DSQRT( PV) 

TF(DARS@@Iy GE, 7, 0)Goule 12 

PHI 1=( DEXP( -(T1**2)/2. ))/SQRT(2. #3. 14159265) 
CALL MDNORD(T1,CDD1) 

DD1=1. 0-CDD1 


54 


12 


13 


14 


15 


20 


3° 


rid 


GO TO 13 

IF(T1. GT. 7. 0)DD1=0. 0 

Pee bl. (=7. 0) )DD1=1.0 

PHT1=0.0 

T2=(X2-Z) /DSQRT( PV) 
IF(DABS(T2).GT.7.0)GO TO 14 

PHI 2=( DEXP( -(T2**2)/2. ))/SQRT(2. *3. 14159265) 
CALL MDNORD(T2 ,CDD2) 

DD2=1. 0-CDD2 

GO TO 15 

IF(T2. GT. 7. 0)DD2=0. 0 

IF(T2. LT. (-7. 0))DD2=1. 0 
PHI2=0. 0 

DLPHA1=( PHI1-T1*DD1)*( DSQRT( PV) ) 
DLPHA2=( PHI2-T2*DD2)*( DSQRT( PV) ) 
IF(DLPHA1. LT. 0. 0)DLPHA1=0. 0 

IF( DLPHA2. LT. 0. 0)DLPHA2=0. 0 

OV 1=( DLPHA1 -DLPHA2 )*D/QR 

RETURN 

END 


*VARIABLES NOT DEFINED IN PREVIOUS APPENDICES 


*VARIABLE NAME DATA 


*MD 
*NAME 
*QM 
*QRR 
*NT 
oe 
zNN 
*NPO 
*TOV 
*QMR 
*SLD 
*OV 
*OV1 
ve 


1 
UICP CONSUMABLES MODEL 

(NOT USED) 

(NOT USED) 

1500 (LIMIT ON THE NUMBER OF ITEMS THE PROGRAM WILL 
PROCESS 

2 (NUMBER OF MEASURES OF PERFORMANCE - SMA & MSRT) 
1 OR O - PRINTOUT OPTION (1=SUMMARY, O=FULL PRINTOUT) 
COMPUTED SMA AND MSRT 

(NOT USED) 

DAYS OF SAFETY LEVEL FOR THE ENTIRE FOUR DIGIT COG 
MSRT FOR EACH NIIN 

EXPECTED NUMBER OF BACKORDERS 
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C --- 
C --- 


10 
ial 


de 


iS 


@ aoc 


APPENDIX H. SUBROUTINE CDFP 


SUBROUTINE CDFP(2ZZ,K,P,C,D) 

ROUTINE TO CALCULATE POISSON CDF AND MASS 
CUMULATIVE POISSON DISTRIBUTION 

REA“ 60222 PP CG CCl DDG 


REAL ZZ,P,D 
INTEGER K,I 

1F( RUT 0)GO 10 12 
ZZZ=ZZ 

PP=DEXP( -ZZZ) 
CC=PP 

Cc i-ona 


IF(K. EO) 0) GO TO 11 
KK=5*IFIX(ZZ+0. 5) 
IF(2Z.GT, 10.0, AND. K. GE. KK)GO ZO 15 
DO 10 I=1,K 
CC1=CC 
PP=PP*Z2Z/DBLE(1) 
CC=CC+PP 

CONTINUE 

P=Pp 

c=CC 

DD=1. 0-CC1 

D=DD 

RETURN 

P=0. 0 

C=0.0 

D=1.0 

RETURN 

P=0. 0 

c=, 0 

D=0. 0 

RETURN 

END 


SUBROUTINE CDFB(ZZ,K,PVAR,C,NB) 
ROUTINE TO CALCULATE NEGBINOMIAL CDF 
REAL ZZ,C,PVAR 

REAL*8 ZZZ,PP,CC,BR,R,BK,S22,B,BQ 
INTEGER K,I,NB 

NB=0 

Z22=22 

$22=PVAR 

BR=ZZZ/S22 

BQ=S22/Z2ZZ 

IF(BQ. LE. 1.0)GO TO 8 

R = 1.0-BR 

BK=( ZZ2**2) /(S22-Z2Z) 
IF(BK*DLOG( BQ). GT. 9.0)GO TO 8 
PP=BR***BK 

CC=PP 

TFCK 2 2G 0) GO TO 17 
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GO TO 9 
8 NB=1 
RETURN 
Geo: I=1.K 
B=DBLE(I-1) 
PP=PP*R*(B+BK) /DBLE( 1) 
CC=CC+PP 
10 CONTINUE 
11 C=CC 
RETURN 
END 
* 
*ALL VARIABLES EXCEPT OUTPUT VARIABLES (P, NB) PREVIOUSLY DEFINED 
bd 
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APPENDIX I. SAMPLE CHURN 


HOA | COG 4a | COG aia yaaa 
| CHURN TPREO! CHURN | EREO Ghar SS a 
eS CE iii. (2020 


‘ 
asf SA 
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a|—~ 


nae 2 ae 
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14. SAMPLE CHUR ONT TED 
mega sat Oe ina COC Tar 1 COG IhuD 
ieee Cie NP REO] CHURN | FREO] CHURN TFREOQ! CHURN | 
a -36 ee ee Eee eee 
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